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INTRODUCTION 


Le problème du plan d'expérience et de son analyse statistique a 
été introduit et étudié par Sir R. A. Fisher. Il a proposé plusieurs ty- 
pes de plans d'expérience ; F,. Yates en a découvert d'autres, en par- 
ticulier le ‘Bloc Incomplet Equilibré'', R.C. Bose et l'école indienne 
ont continué l'étude des blocs incomplets et ont introduit le ‘'Bloc In- 
complet Partiellement Equilibré"" qui sera l'objet de cette étude. 


Dans le domaine des recherches agricoles, de la biologie ani- 
male, les plans d'expérience sont d'un emploi courant : pour comparer, 
par exemple, le rendement de diverses variétés de blé plantées dans 
un terrain connu, ou l'effet de divers traitements sur la quantité de lait 
produite par des vaches appartenant à plusieurs races, ou encore l'ef- 
fet de divers régimes alimentaires sur de jeunes bestiaux, etc. 


Nous emploierons la terminologie de R. C. Bose qui est celle 
empruntée aux problèmes de recherches agronomes : 


Le terme variété désignera les quantités dont on expérimente les 
effets, variétés de blé ou traitements alimentaires dans les exemples 
donnés. 


Le terme bloc désignera le groupe d'unités sur lequel agissent 
les variétés, bloc de terrain divisé en parcelles ou groupe d'animaux. 


Un des premiers plans d'expérience introduits est le "Bloc com- 
plet'' oùutoutes les variétés sont réparties au hasard dans chaque bloc. 
Ce plan d'expérience nécessite dans un bloc autant de parcelles que de 
variétés, donc un grand nombre d'expériences et, quand le nombre de 
variétés est élevé, un grand nombre de parcelles dans chaque bloc ce 
qui est parfois impossible surtout dans le domaine de la biologie ani- 
male. 


De plus il n'existe pas d'analyse mathématique rigoureuse de ce 
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plan d'expérience (il n'est pas tenu compte par exemple de la variation 
de fertilité du sol). 


Pour toutes ces raisons d'autres plans d'expérience lui sont 
substitués la plupart du temps. 


Les carrés latins orthogonaux permettent de construire des plans 
d'expérience préférables aux précédents. 


Leur analyse est rigoureuse et permet de tenir compte de la va- 
riation de qualité des blocs et parcelles. 


Cependant il est parfois impossible de l'utiliser puisque le pro- 
blème de la construction des carrés latins orthogonaux n'est pas résolu 
pour tous les nombres. 


Chaque variété est répétée autant de fois qu'il y a de variétés à 
expérimenter, ce qui nécessite beaucoup de temps et de dépenses quand 
le nombre total de variétés est élevé. 


Très souvent, dans les expériences biologiques en particulier, 
le nombre de parcelles de chaque bloc est imposé. C'est le nombre 
d'animaux d'une même portée par exemple. 


I1 est donc nécessaire d'utiliser un plan d'expérience où chaque 
bloc ne contient pas toutes les variétés. Les ‘'Blocs Incomplets'' ont 
été introduits pour remplir ces conditions. 


Tout d'abord le "Bloc incomplet équilibré" qui est le plus simple 
et le plus pratique : l'analyse en est facile et la comparaison des deux 
variétés a la même précision pour n'importe quel couple, le nombre 
de répétition de chaque variété est laissé libre. 


Cependant la construction de plan d'expérience de ce type n'est 
pas donné par un procédé formel : il n'existe aucune règle exhaustive 
de construction. Nous ne connaissons pas de groupe de conditions né- 
cessaires et suffisantes pour construire un bloc incomplet équilibré. 


De plus l'emploi du bloc incomplet équilibré est parfois incom- 
mode : Sir R. À. Fisher a démontré (Ann. of Eugénics 1940) que la 
construction d'un bloc incomplet équilibré n'était possible que si le 
nombre de oloc était au moins égal au nombre de variétés, ce qui né- 
cessite un grand nombre d'expériences. De plus, pour certaines va- 
leurs du nombre de variétés, la construction d'un bloc incomplet équi- 
libré exige un trop grand nombre de répétitions : par exemple, pour 
huit variétés et des blocs de trois parcelles, le plus petit nombre 
de répétition est vingt et un. 


Te 


qe 
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Pour lever ces difficultés, R. C. Bose et Nair introduisent en 
1939 une nouvelle catégorie de plan d'expérience : le bloc incomplet 
partiellement équilibré à deux classes associées, qui comprend comme 


cas particulier le bloc équilibré. En 1952, Bose et Shimamoto donnent 


la définition la plus générale, du bloc incomplet partiellement équili- 
bré à m classes associées, 


L'objet essentiel de cette thèse est d'étudier des constructions, 
des conditions d'existence de tels plans d'expérience et d'établir l'ana- 
lyse de la variance dans le cas le plus général, 


Dans un chapitre préliminaire nous rappellerons, avec la défini- 
tion, toutes les relations existant, dans le cas le plus général, entre 
les paramètres du plan d'expérience. 


Le premier chapitre sera relatif à la construction dans le cas où 
le nombre de classes associées est supérieur à deux. Nous générali- 
serons une construction donnée dans le cas de deux classes associées, 
ayant pour base les géométries projectives finies. Nous l'étudierons 
endétaildans le cas de la géométrie à trois dimensions construite sur 
le corps de Galois à 4 éléments. Ensuite nous donnerons une méthode 
de construction de plan d'expérience cyclique et de plan d'expérience 
quasifactoriel à n dimensions. Û 


Les deux chapitres suivants étudieront le cas particulier de deux 
classes associées. 


Dans le chapitre deux, nous rappellerons les travaux de Connor 
et Clatworthy sur la matrice d'incidence d'un bloc incomplet partielle- 
ment équilibré, puis nous établirons des conditions d'existence de tels 
plans d'expérience dans le cas où le nombre de blocs est inférieur au 
nombre de variétés, 


Dans le chapitre trois, nous dresserons la liste de tous les blocs 
présentant quatre répétitions ; certains figurent dans la classification 
établie par Bose Clatworthy et Shimamoto, d'autres étaient encore in- 
connus. Pour ces derniers nous donnerons les schémas d'association, 
il nous est apparu au cours de ce travail qu'il pouvait exister d'autres 
schémas que ceux envisagés par les auteurs de la classification ; nous 
les décrirons. 


Dans le chapitre quatre, nous établirons suivant la méthode gé- 
nérale donnée par Rao, l'analyse de la variance "intra-inter bloc" et 
la valeur du coefficient d'efficacité, en mettant en évidence les fonc- 
tions des paramètres de seconde espèce qui interviennent dans cette 
analyse. Puis, sous certaines conditions, nous établirons l'analyse 
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intra bloc" dans le cas où tout un bloc manque dans le plan d'expé- 
rience, ce cas pouvant se présenter dans la pratique quand le bloc est 
une unité expérimentale. 


Les principaux résultats ont été donnés dans 3 notes aux Comptes 
Rendus de l'Académie des Sciences. 


Je suis heureuse de pouvoir exprimer ici toute ma reconnaissance 
à M. Dugué sans qui ce travail n'aurait jamais vu le jour ; pour les 
très nombreux entretiens qu'il a bien voulu m'accorder, pour ses nom- 
breux conseils, ettoutes les améliorations qu'il a apportées à ces tra- 
vaux. 


Que M. G. Darmois, dont les encouragements m'ont été siutiles 
et qui a bien voulu me faire l'honneur de présider ce jury de thèse, 
veuille trouver ici l'expression de toute ma gratitude. 


Je remercie également M. Garnier qui a accepté de me proposer 
un sujet de seconde thèse et de faire partie du jury. 


CHAPITRE PRÉLIMINAIRE 


Pour comparer v variétés dans des plans d'expérience formés de 
blocs de k parcelles, où k est inférieur à v, nous disposons de deux 
types de plans d'expérience : 


- les “Blocs Incomplets Equilibrés" ; 
- les "Blocs Incomplets Partiellement Equilibrés". 


Le Bloc Incomplet Equilibré est le plus simple des deux et le plus 
efficace lorsqu'il existe. Nous allons donner un exemple de plan de ce 


type. 


Considérons 6 variétés, 1, 2, 3, 4, 5, 6 réparties en blocs de 3 
parcelles de la façon suivante : 


BB BB O1 ©  « 
D OO D où 


2 
2 
2 
3 
3 


M OM ke M 
OO & © NN NN 
D OO où À «wo 


Nous remarquons que chaque variété apparaît 5 fois, et que cha- 
que couple de variétés apparaît dans 2 blocs. 


Les plans d'expérience possédant ces propriétés : 


a) le nombre de répétitions de chaque variété est indépen- 
dant de la variété, 


b) le nombre À:; de blocs contenant le couple de variétés 
vi v; est indépendant de i et j, 


sont appelés ‘Blocs Incomplets Equilibrés". 
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Nous allons donner ladéfinition générale en utilisant les notations 
suivantes : 


y nombre de variétés, 


b 1 "n blocs, 

k ul "1 parcelles dans chaque bloc, 

r ul 1 répétition de chaque variété, 

À L "blocs contenant un couple de variétés. 


Définition d'un Bloc Incomplet Equilibré. 


Un bloc Incomplet Equilibré est un arrangement de v variétés en 
b blocs de k parcelles chacun, tel que : 


1/ Chaque variété apparaisse au plus une fois dans chaque bloc. 
2/ Chaque variété est répétée r fois. 


3/ Chaque variété v; apparaît avec la variété v; dans À blocs, 
étant indépendant de v; et v;. 


Relations entre les paramètres. 


Les paramètres v, b, r, k, À, ne sont pas indépendants. 


Le nombre d'expériences faites est égal au nombre total de par- 
celles : bk, ou au nombre total de variétés distinctes ou non : vr. 


- Chaque variété v; apparaît dans r blocs. Ces r blocs con- 
tiennent r(k-1) variétés différentes de v., de plus v; apparaît avec cha- 


cune À fois : 
r(k-1) = À (v-1) 


- Sir R. À. Fisher a démontré qu'il était nécessaire que 
b>v 


et R.C. Bose, dans le cas où b est de la forme b = rn et où le bloc est 
résolvable", a démontré l'inégalité plus forte : 


b>v+r-1 
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Ces conditions ne sont pas suffisantes pour assurer l'existence 
d'un telplan d'expérience. Un groupe de conditions nécessaires et suf- 
fisantes reste encore inconnu. 


Nous allons rappeler les principales méthodes de constructions 
de Blocs Incomplets Equilibrés. 


Constructions de Blocs Incomplets Equilibrés. 


1/ Unpremier groupe de Blocs Incomplets Equilibrés est obtenu 
en prenant comme variétés les points d'une géométrie projective finie 
de dimension m construite sur le corps de Galois à p" éléments (p nom- 
bre premier), et comme blocs les géométries projectives finies de 
dimensions, s < m, contenues dans la précédente. 


Il est facile de démontrer que les paramètres ainsi définis véri- 
fient les relations et conditions précédentes, 


2/ Un deuxième groupe de Blocs Incomplets Equilibrés est obte- 


nu, de façon analogue, à partir des géométries euclidiennes finies. 

3/ Bose a donné une autre construction, utilisant les groupes 
abéliens additifs ; il a pu ainsi établir une liste systématique et presque 
exhaustive de ces plans d'expérience. 


Les Blocs Incomplets Partiellement Equilibrés sont une généra- 


lisation des Blocs Incomplets Equilibrés. 


Ils existent dans des cas où il n'existe pas de Blocs Incomplets 
Equilibrés, par exemple pour b < v ; ou bien ils existent avec un nom- 
bre de répétition inférieure à celui du bloc équilibré, par exemple pour 
v=8, k=3, la plus petite valeur de r est 21 dans le bloc équilibré, et 
3 dans un bloc partiellement équilibré. 


R. C. BoseetK.R. Naïir ont donné les premiers la définition d'un 
Bloc Incomplet Partiellement Equilibré à 2 classes associées seule- 
ment. Nous allons l'énoncer et l'illustrer d'un exemple. 


Définition. 


Un Bloc Incomplet Partiellement Equilibré à 2 classes associées 
est un arrangement de v variétés en b blocs de k parcelles chacun, tel 


que : 
1/ Chaque variété apparaisse une fois au plus dans chaque bloc. 


2/ Chaque variété apparaisse exactement r fois dans l'ensemble 
du plan d'expérience. 
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3/ I1 existe une relation d'association.entre 2 variétés quelcon- 
ques, satisfaisant aux conditions suivantes : 


a) Deux variétés v; et v, sont premières ou deuxièmes as- 
sociées. 

b) Chaque variété a n, premières associées et n, deuxiè- 
mes associées; n, et n, étant indépendants de la variété 
initiale. 

c) Deux variétés ièmes associées apparaissent ensemble 
dans À;blocs (i = 1,2). 

d) Soient deux variétés ièmes associées le nombre de va- 
riétés ièMeS associées de la première et jèmes asso- 


ciées de la seconde est Pi; (u,"3, 7 271, 2). Ce nombre 
est indépendant du couple initial, 


Ilest évident que nous avons là, une généralisation du bloc équi- 
libré. En effet, quand À, = À, les deux classes d'association se con- 
fondent, et nous retrouvons la définition du bloc équilibré. 

Relations entre les paramètres. 


En raisonnant comme pour le bloc équilibré, nous obtenons les 
relations suivantes : 


ven; tnt! 


vr = bk 


n,À,+n, À,= r(k-1) 


De plus, n, étant une constante 


TD 
(l] 
T 


Exemple. 


Considérons le bloc incomplet partiellement équilibré, à deux 
classes associées, ayant pour paramètres : 


Vo tOu=ND kr 
Die 60 O0 EN RAA ENS EN = 0 
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et pour paramètres p 
1 = 1 = 1 = = = 2 = 
PR D OP HI UP. 0,4p *"1 


le tableau d'association est le suivant : 


variété premières associées deuxièmes associées 
1 2 3 4 SRG 7 89 
8 2 3 4 5.6 a 1 9 
9 2 3 4 GS 6 7 1 8 
2 1 3 4 5 CIO 6 7 
6 1 3 4 5 OO 2 7 
% il 3 4 5 & À 7 € 
3 1 26 % # 4 5 
4 il 2 G 7 # € ë 5 
5 1 2 RG Ü 8 9 3:11a 4 
Le plan d'expérience a pour schéma : 

1 2 3 6 SES fl 9 $ 

1 6 4 6 9 5 2 8 5 

1 % 5 fl 8 4 PAM 


Nous pouvons vérifier qu'il possède les propriétés énoncées dans 
la définition. 


R.C. Bose et K.R. Naïir ont généralisés ces plans d'expérience 
en considérant un nombre quelconque de classes d'association, 


Définition générale d'un Bloc Incomplet Partiellement Equilibré. 


Un Bloc Incomplet Partiellement Equilibré à m classes associées, 
est un arrangement de v variétés en b blocs de k parcelles chacun, tel 
que : 


1/ Chaque variété apparaisse au plus une fois dans chaque bloc. 


2/ Chaque variété apparaisse exactement r fois dans l'ensemble 
du plan d'expérience, 


3/ 11 existe une relation d'association entre 2 variétés quelcon- 
ques, satisfaisant aux conditions suivantes : 
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a) Deux variétés v; et v; sont premières, deuxièmes, ... 
ou mièmes associées. 


b) Chaque variété à n, uièmes associées, n, étant indépen- 
dant de la variété initiale, 


c) Deux variétés ièmes associées apparaissent ensemble 
dans i blocs. 


d) Soient deux variétés uèmes associées, le nombre de va- 
riétés ièmeS associées de la première et jèmeS asso- 
ciées de la seconde, est P', 


GEMEMNE GR PE don das) 
Ce nombre est indépendant du couple initial. 


Les nombres v, b, r, k, n;, À;, sont appelés paramètres de 
première espèce, ce sont des entiers positifs. 


Les nombres pi, paramètres de deuxième espèce, ce sont des 


entiers positifs ou nuls. 


Relations entre les paramètres, 


Le nombre d'expériences faites est égal au nombre total de par- 
celles : bk, ou au nombre total de variétés distinctes ou non : vr 


- Le nombre de variétés apparaissant dans les mêmes 
blocs que v, peut être calculé de deux façons différentes : 


Et 
- n, étant une constante : 
Pi, x P'; (i, 15 Us 4, 2, 10) 
Nu Pi; = np}; 5 np; 


- Le calcul dunombre de variétés uèmes associées fournit 
la relation : 
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rées, 


= 
| Dr Dé ss De (ls 
ij ij 1J 


M | 


Em 
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HE:0 siu pi 


SIEUR 21 


Remarquons que l'ensemble des paramètres de seconde espèce 
p;. peuvent être écrits sous la forme de m matrices symétriques car- 


.. » (M) 


Ces notations sont celles de R.C. Bose, elles sont employées 
dans tous les articles ayant trait à ces plans d'expérience. 


s |] 


CHAPITRE | 


CONSTRUCTION DE PLANS D'EXPERIENCE A PLUS 
DE 2 CLASSES ASSOCIÉES 


I - CONSTRUCTION UTILISANT LES GEOMETRIES PROJECTIVES 
FINIES - 


Bose et Nair ont donné, en 1939, une construction de blocs in- 
complets partiellement équilibrés à deux classes associées utilisant 
comme variétés les points d'une géométrie projective construite sur 
le corps de Galois à p" éléments (p premier). Cette construction est 
elle-même une généralisation de la construction de blocs incomplets 
équilibrés. 


Nous allons donner ici, une construction dans le cas où le nombre 


de classes associées est 3, et où la dimension de la géométrie pro- 
jective utilisée est 3, 


1/ - Construction générale. 


Considérons le corps de Galois à p" éléments. Un point d'une 
géométrie projective finie à m dimensions GP(m, p") est un ensemble 
ordonné de m+ 1 éléments du corps de Galois non tous nuls. 


Deux ensembles représentent le même point (8,59. 106888) 
(60-008 DOS OS, ON Me umo 


b #0 ue (p') 


La droite qui joint deux points p, et p, est définie comme l'ensem- 
ble des points Hp, + H,Pp, 


HE (pu,#0 


Une GP(m,p')a 


IPESDIEER TP RE RDREDOInLS: 
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Chaque GP(m, p") contient exactement : 


CPP EEE D) ET (DOUÉ. 21e 
EDS D D). "D: 


GP(Kk, p') 


Chaque GP(k, p°) d'une GP(m, p") donnée est contenue dans 


(PEARL ES pr) PRET. + pr) 


GP(s, p° 
(pie L alt se p°') Ve p” ( P ) 


avec s > k 


Si k = 0 nous obtenons ainsi le nombre de points d'une GP(s, p') 
contenus dans une GP(s, p") donnée, 


Si k = 1 le nombre de droites d'une GP(m, p") contenues dans une 
GP(s, p'). 


2/ - Construction de blocs incomplets partiellement équilibrés à 


3 classes associées, 
Considérons une géométrie projective de dimension 3 construite 
sur le corps de Galois à p" = a éléments, GP(3, a). 


Faisons choix d'un plan |], et dans ce plan, faisons choix d'une 
droite D et d'un point À non situé sur D. Remarquons que n'importe 
quel point de Îl, non situé sur D, peut jouer le rôle du point A cela 
revient à faire une permutation sur les variétés. Nous considérerons 
les points de GP(3, a) non situés dans [|] comme des variétés 


v=l+a+a?+a-(lrara) = a 
Nous identifierons les blocs avec les plans qui ne passent pas par D. 


milratasna) (Grara) (are) 


- (1+a) = a? + a? 
(1+a+a°?) (a+a°) a? 


Chaque bloc 1 + a + a points moins les points du bloc situés 
dans N 


ke: litia tra (1h a), =187 


Chaque point est situé dans r plans dont un passant par D 
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(a+a?+ ai) (a?+ a?) : : 


(a+ a?) a? 


=za+a 


Ces points forment un plan d'expérience, en effet vr = bk. 


Définition des classes d'association. 


Un point M sera premier associé des points de la droite MA 


heal 


Il sera deuxième associé des points du plan MD 


Troisième associé des autres points du plan d'expérience 


= 


Deux points premiers associés apparaissent ensemble dans 


Al /r anpiocs 


Deux points deuxièmes associés apparaissent ensemble dans 
À, = a blocs 
Deux points troisièmes associés dans 


A3 =1+a blocs 


Le plan d'expérience ainsi formé est un bloc incomplet partiel- 
lement équilibré à 3 classes associées : 


NRA STE À,n, = r(k-1) 


Les paramètres de seconde espèce sont donnés par les matrices: 


ee eo, KO no 


OMS ON 
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0 0 0 
0 n P;. = 0 
n, n,-n, n. 
O1 n, 
p} = 1, 0 n, 
Nn;, N, n,-n,-n,-1 


3/ - Exemples de construction. 


Cette géométrie GP(3, 2) est formée de 


(es NN 


(0) 
0 
1 
0 


(0) 
0 
(9) 
1 


- À sera le point p , et D la 


put LR OMC Pa OS 
p,<1 OMEIMO pre © 0 
P,= 1 OO Pons is 
sa 1OSTALO Pan 


droite P, P, 


- les points de [1 a exclure sont : PP PP, 


207 


OR TERT 
Lee 
LIL 
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- le point premier associé de p est p+p, 


- les 3 points seconds associés de p sont : 


PRE SP RER 
PRSD ET 
DSC 


Les droites de cette GP(3,2) sont : (nous remplacerons le point p 
par i) 


LL 2 2. & 6 sin Gin AIT GANT 12: 
SG 2" 4 9 3) "7 215 SCIE CMS 
lee? nn (ES LI SM OT T Hs a € 70 TOI T 
1 6 aÙl in an 9201220015 5 10 15 OO T0 
bn Re 7 2 10 14 7 5 13 14 SNS 
ll jo) 16 2 13 15 4 6 13 Ge NT JO CIEL 5 LE 
1 14 15 SEL O 4 88 14 COR 15 JON PSS? 


Nous construirons les plans de la géométrie en prenant les cou- 
ples de droites passant par un point i. Nous ne conserverons de ces 
plans que ceux qui ne passent pas par la droite 1, 2 et nous enlèverons 
lesipointside NPASOten SRG RS 2 5 


4 7 T2 


LOIS AS 
4 TON 
CPS GRIS 465 a 
4 Cnil 
4 CL OS AE 
6 A0, A2 5 
LS 9 14 
4 12 13 14 
{l SOLS 
CAT OST 


Tableau d'association. 


- Le point premier associé du point p est le point p + pe 
c'est-à-dire p + 100 0, 
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PMAOMOMOETMBA pour premier associé A0: "0 = P, 

Ps OT OISE = - TROUS P,, 

DR OMIS = - - CL OIE RE P., 

DU lle ser = = LT P,, 

- Les trois points deuxièmes associés de p sont : 

EE CAES LA on 2 prtq+t 
pour p, ce sont ORLEOMIE P, » OROMTSITE= P,,» D TRTETES ps 
pour p, ce sont LSLROPIE pe ADM LEE P3 LOL SUIS P, 


Nous pouvons résumer ces résultats en un tableau rectangulaire 
où chaque variété i sera : 


- première associée de la variété située sur la même 
rangée que i; 

- deuxième associée des trois variétés situées sur la mé- 
me colonne que ji; 


- troisième associée des variétés non situées sur la même 
rangée ou la même colonne que i. 


Remarque. 


De façon générale le nombre des variétés étant a’, chaque va- 
riété i ayant a - 1 premières associées et a? - 1 deuxièmes associées, 
le tableau d'association des variétés se présentera sous forme d'un 
rectangle ayant a variétés dans une rangée et a? variétés dans une co- 
lonne. Les variétés d'une même rangée sont premières associées cel- 
les d'une même colonne deuxième associées. 


II - CONSTRUCTION A L'AIDE DE TREILLIS CIRCULAIRES - 
A - CONSTRUCTION DES TREILLIS - 


Dans un espace vectoriel à s dimensions sur les corps des nom- 
bres réels R. choisissons un sous-espace de dimension s-2 P,, 
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et n sous-espaces de dimension s - 1 ayant en commun P,, 
Pair et MR. =P,, 
Continuons la construction de cette chaîne : 
- dans P., faisons choix d'un sous-espace P,., 
ps CE. 


- dans chaque R.., faisons choix de n espaces de dimension 
s - 2 contenant pi: 


R:> (æ R: CZ R, 


Finalement choisissons une droite p}; 


PPT ec APE CR: 


et dans chaque R, obtenu par le processus précédent, choisissons n 
plans R, passant par la droite P, 


ReNG Ci 
MR, =p, 


Dans chaque plan R, construisons un treillis circulaire d'ordre 
n, c'est-à-dire n cercles concentriques et n diamètres. 


Le treillis ainsi formé est composé de 2 points, chaque cercle 
et chaque diamètre, renfermant 2n points. 


Si nous considérons les cercles et les diamètres comme des 
blocs, les points comme des variétés, nous obtenons un plan d'expé- 
rience avec les paramètres suivants : 


Montrons maintenant qu'il s'agit 1à d'un bloc incomplet partiel- 
lement équilibré à s+1 classes associées, 
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B - DEFINITION DES CLASSES D'ASSOCIATION - 


Choisissons une variété v les autres se répartissent en s +1 
groupes : 


1/ les variétés apparaissent avec v à la fois sur le même dia- 
mètre et sur le même cercle : ilyena n,= 1. 


2/ les variétés apparaissant soit sur le même cercle soit sur le 
même diamètre que v, et n'appartenant pas à la première catégorie ; 


ilyenan,=4{(n-1). 


3/ les variétés appartenant au même plan que v (sans appartenir 
aux autres classes) 


n, = 2n°-4n+2 


n, = 2(n-1) 
4/ les variétés apparaissant dans le même R;, que v 


ni ?n"2n? 


D 
(I 


2n°(n-1) 


S+ 1) les variétés de R,, n'appartenant pas au même R.,que v 


n D Pine Se © Pie, Ar 


S+1 
nn = 2n°i (n-1). 


Conformément aux notations de R.C. Bose, appelons À le nom- 
bre de blocs dans lesquels deux variétés ièMES associées apparaissent 
ensemble, 


i=1 


r(k-1) = n'_À, = 2(2n-1) 


i=] 
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La loi d'association est telle que les paramètres de seconde es- 
pèce p°. sont bien indépendants du couple initial choisi. Les matrices 
des paramètres de seconde espèce s'écrivent : 


AN 
ue 
=) 
1 
CEA 
D 4 
e) 
©Q [æ)] © 
—) SOS 


ses ee ee © © + e 


0 
2(n-2)  2(n-1) 0 
0 


T 
(l 
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Valeurs qui vérifient le second système de conditions : 


Ve Di ASTON =) 
SD E DONS 
Cat MIO Er NS 


Nous avons ainsi construit un bloc incomplet partiellement équi- 
libré à s+1 classes associées, 


CEICAS PARTICULIERS. 


Nous allons donner les paramètres de plans d'expérience qui 
peuvent être d'un usage courant. 


[==2] C. Radhakrihna Rao en a donné des exemples dans la revue 
Sankya. 


Exemples de construction, 


Dans chaque treillis les variétés sont numérotées sur les diamè- 
tres à partir du cercle intérieur et en tournant dans le sens de rotation 
positif, 

1/-8=3, n=2,. 

Nous écrirons, sur le schéma, les numéros des variétés du pre- 
mier plan droits, les numéros des variétés du second plan seront en 
italique. 
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Diamètres Cercles 
(I, 3,9) premier plan CPS MST) 
(Es Opus (2, 4,6, 8) 


CSAIL, 13191 
(10, 12, 14, 16) 


(9/10: -11,012) 
(13, 14, 15, 16) 


second plan | 

Les couples premiers associés sont : 

(5-3) (24) (6527)(6;%0) (910) (10,22) (13° 15) (245716): 

Les variétés secondes associées de la variété 1 sont : 
2,45, 

Les variétés troisièmes associées de la variété 1 sont : 

678; 

La simple lecture du treillis fournit les variétés associées. 

2/-s8=3, n =3, 

25419, 47:59 000) 


En) AE ei 12) | premier plan | 
145215; 167-177 18) 


14 27, 10518: 
y— 058, 11,14; 
6 PhONeLA FES; 


nn ne fn 
© ND 


20, 21, 22, 23, 24) (19/%22,225,,28 31 
26,227, 20/29/2350) second plan 20 23, 26, 29, 32, 
32, 33, 34, 35, 36) (21, 24, 27,30, 33, 


38, 39, 40, 41, 42) (37, 40, 43, 46, 49, 
44, 45, 46, 47, 48) troisième plan je 41, 44, 47, 50, 
50, 51, 52, 53, 54) (39, 42, 45, 48, 51, 


16) 
17) 
18) 


34) 
35) 
36) 


52) 
53) 
54) 
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18 26 


20128 
3/-s=4, n = 2, 
Diamètres Cercles 

OS En) premier plan (ALES OR T) 
(5 #6, 7, 0) (ERA AGE) 

ler R, 
(19240, 411,:,12} } deuxième ns OH Ris es 
(11400 rS 216) (0412 44016) 
HIT 018 10 20) ds (172%199721,%23) 
(21952223 24) Pl P (18, 20, 22, 24) 

2ème R, 
(25, 26, 27, 28) | j.jame plan Fée CPC NE A) 
(29, 30, 31, 32) P (26, 28, 30, 32) 


4/ - Construction générale : s = 2n. 


Diamètres : (names tn (ri) 


ler plan 1-0 0, nel 
LEP sata rersserssnaesscocscsesds se 
avec : premier R, nème plan: dent) #22 nl 
ler plan i= (n-1)n°,..., (n-1)(n°+1) 
OS Re un one se une no dise VC 


nème plan : 1e une 
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Cercles : (ARS CGT EME) 
dans le premier plan : JE dh) 660 0 
dans le tè plan: j= 2{(t-l)n°+1,... 2 bn mn 


II - CONSTRUCTIONS DE PLANS D'EXPERIENCE QUASIFACTO- 
RIELS DE DIMENSION n - 


Un bloc incomplet partiellement équilibré est dit quasifactoriel 
de dimension n, lorsque le nombre de variétés v est de la forme 


VS pa D: Che D) 


1 2 


DEND cp: étant des nombres quelconques. Les plans d'expérience 
quasifactoriels de dimension deux, ont été étudiés par Radhakrishna 
Rao. Nous allons faire ici une étude générale qui permettra d'augmen- 
ter le nombre de classes d'association. 

A - DEFINITION DES CLASSES D'ASSOCIATION - 


Dans un espace vectoriel à n dimensions ayant pour base l'en- 
semble des vecteurs linéairement indépendant 


les v variétés peuvent être rangées en réseau formant un hyperparal- 
lèlépipède aux côtés parallèles aux vecteurs e, : ily a p. variétés sur 
chaque parallèle au vecteur e.. 


1/ Deux variétés seront i° associées, (i = 1, ..., n) si elles se 
trouvent sur la même parallèle au vecteur e:;, d'où 


nas al 


2/ Deux variétés seront j° associées, si l'espace à 2 dimensions 
construit parallèlement aux#vecteursres ete (,-:k= 1 nm) 1er 
passant par l'une, contient l'autre, (sans qu'elles soient i ou k asso- 


ciées) 
nj = (p;- 1) (px - 1) 


Il y a Ci classes d'association de cette catégorie. 


3/ Deux variétés seront u° associées si l'espace à 3 dimensions 
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construit parallèlement aux vecteurs e;, e;, e,, contenant l'une con- 
tient l'autre (sans qu'elles appartiennent aux deux premières catégo- 
ries), 


ne (Di l}(pie st) (pri) 


Je = et costs 2 LE 
iCjCk 


Ily a (e classes d'association de cette catégorie. 


Les classes d'association construites sont bien disjointes puisque 
nous avons considéré chaque ensemble de variétés comme ouvert. Il 
y a en tout 


NC HC He. 6. 


n 


N=12 29 1 classes d'association. 


B - BLOCS INCOMPLETS PARTIELLEMENT EQUILIBRES A 
N CLASSES ASSOCIEES - 


Supposons qu'il existe n blocs incomplets Equilibrés ayant pour 
paramètres : 


VD TT, k- À HA (Eee Mn) 
(ils ont le même nombre de parcelles). 


Considérons le réseau des variétés comme forme d'ensembles 
de p, variétés disposées parallèlement au vecteur e;, et construisons 
pour chaque ensemble un bloc incomplet équilibré de paramètres p,, b,, 
hr ANR AU 


En juxtaposant les n groupes de ml P; (j # i) blocs incomplets 
E 


équilibrés ainsi construits, nous obtenons un plan d'expérience dont 
les paramètres sont : 
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(; -1)(p; - 1)...(p,- 1) 


. , N) avec N=2" - 1 


Ce plan d'expérience est un bloc incomplet partiellement équili- 


bré à N classes associées, en effet : 


n n 


k » (b; n P;) F (> r:) n P; Car kb; SAC: D: 
i= j#i i=1 
EE - 
N n 
SEN ANS Te (k-IJENLr,, Car AID -1) = Fr (KP) 
S=] Î=1 


Lorsque les p, sont distincts, les seules valeurs pratiquement 
intéressantes de n sont 2 et 5. 


Remarquons cependant, que chaque fois que deux ou plusieurs p, 
sont égaux il y a réduction du nombre N de classes d'association :n 
peut alors prendre des valeurs supérieures à 3. 


CICASLPARTICULIERSS 


1/ Nous pouvons toujours prendre comme plan d'expérience équi- 
libré le plan de paramètres : 
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Paramètres du plan Paramètres du plan 
équilibré quasifactoriel obtenu. 


mans ner 


4(1+3°) 
6+10. 3° 
7(1+3°) 
9(1+4.3") 
4(1+4°7) 
7(1+4°) 
8(14. 4" 
4(1+4°) 
5(1+4°) 
53.5" 


3 
5 
3 
4 
2 
4 
7 
3 
4 
3 


2 
2 
1 
il 
2 
2 
3 
3 
3 
3 


3 
3 
3 
3 
4 
4 
4 
4 
4 
5 


4 
6 
Ü 
9 
4 
7 
8 
4 
5 
5) 


CR TRS OS EP CO COM CO ECS 


En le juxtaposant, suivant la méthode exposée au plan d'expé- 
rience équilibré : 


Vas qd L'HAES se LA À ESS KXAIp 


Nous obtenons un bloc incomplet partiellement équilibré à 2n+1 classes 
d'association : 


PRG De DD Hd CK =D, TM, (A; Azml, te 0 
Nous pouvons calculer ainsi le nombre de classes d'association : 


- pour la première catégorie n classes se confondent, il 
n'en reste que deux distinctes ; 

- il en est de même pour les n-1 autres catégories ; 

- la dernière ne comporte qu'une classe, 


Nous allons dresser la liste des principaux plans d'expérience 
qui peuvent être ainsi construits. 


2/-Cas:n=3, v=p.p.P. 
En rangeant les variétés dans un parallèlépipède, nous obtenons 
un plan quasifactoriel à N = 2? - 1 = 7 classes associées. 
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Lorsque deux des trois plans utilisés.équilibrés sont identiques 
il se produit la même réduction de classes d'association que dans le 
cas précédant, 


Nous donnons la liste de ces principaux plans d'expérience : 


V=P,-P,-P,, b=b,p,p,+b,p,P,+b;p,P,, r=r,+r,+#r,, K, À, IA AS NEO 
He47-- 0) 
Paramètres 


Paramètres des plans équilibrés du plan 
quasifactoriel 


AABBADABAARMANLAAMAARMARRLARIBANRERABRAAAAR 


gl 
# 
$ 
# 
# 
# 
Ü 
ô 
ÿ 
ÿ 
y 
# 
# 
ÿ 
Ü 
$ 
ÿ 
y 
# 
y 
ô 
ÿ 
; 
; 
j 
f 
# 
# 
# 
Ü 
5 
Q 


3/ - Calcul de la réduction du nombre de classes d'association 
lorsque plusieurs classes se confondent en une seule. 


Supposons que le nombre v puisse se mettre sous la forme 
VS Cl coco avec s+uz=n 
Dans la première catégorie il y a 1+ ( classes. 
Supposons s > u, dans la seconde il y en a : 


ICE 
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dans la troisième 1+ Ci rC?+ C3 


dans la u° VÉMORTNNIEN SE 


dans les s-u suivantes 1+C'+... + C* 


dans la s + 1° (CAE ARE el 

dans la s + 2° (DE RE et NO 

dans la s+u°=n° C: 
Finalement il y a SSH) AC FC ) 


N = s + (s+1) (2° - 1) 


Supposons maintenant s<u 
dans la première catégorie : 1+C, 


ee 


dans la s° catégorie : 1 + ER pe C2 
dans la s + 1° Cros iC are 
dans la u° Craie + (oi 
dans la n° (Bo 


En tout, comme dans le cas s>u 


N=s+(s+1)(2"-1) 


CHAPITRE II 


CONDITIONS D'EXISTENCE DE BLOCS INCOMPLETS 
PARTIELLEMENT ÉQUILIBRÉS À DEUX CLASSES ASSOCIÉES 


Dans ce chapitre, nous allons étudier le cas particulier, très 
utilisé dans la pratique, où le plan d'expérience ne comporte que deux 
classes associées, 


Nous allons tout d'abord rappeler les résultats obtenus par Connor 
et Clatworthy sur la matrice d'incidence d'un bloc incomplet partielle- 
ment équilibré à deux classes associées, 


Dans le deuxième paragraphe, nous déduirons un groupe de con- 
ditions nécessaires et suffisantes d'existence de tels plans d'expérience 
dans le cas où b < v. (Ce cas est très intéressant pratiquement comme 
nous l'avons expliqué dans l'introduction). 


Nous donnerons finalement une méthode pratique pour dresser la 
liste de tous les blocs incomplets équilibrés existant dans ce cas. 


Dans le troisième paragraphe, enfin, nous résoudrons le systè- 
me d'équations diophantiennes permettant de calculer les paramètres 
de première et de seconde espèce de tous les blocs incomplets partiel- 
lement équilibrés comportant 4 répétitions. Lorsque le plan d'expé- 
rience obtenu a déjà été construit, nous donnerons pour référence le 
numéro du plan dans la Table donnée par Bose R C, Clatworthy et 
Shrikhande, Nous insisterons sur les plans d'expérience encore non 
construits, et sur les cas d'impossibilité. Ce travail a été fait pour 
deux répétitions par Roy et Laha. 


À - ETUDE DE LA MATRICE D'INCIDENCE D'UN BLOC IN- 
COMPLET PARTIELLEMENT EQUILIBRE - 


1/ - Cas général, 


La matrice d'incidence d'un plan d'expériences est : 


INSE (n;;) 
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avec n;,= 1 ou 0 suivant que la ième variété (sel. Mev)iapparait ou 
non dans le jème bloc (j=1,..., b). 


Dans le cas d'un bloc partiellement équilibré, chaque variété est 
b 
répétée r fois, donc à An;; = r (i=1,...,v)et chaque variété ap- 


. J=1 
parait dans À,blocs avec ses sèmes associées, 


b 


ps UE (AIT dde 1; SO) V) 


u=] 


Appelons N' 1a matrice associée de N 


et M la matrice obtenue en remplaçant dans NN' r par la variable z: M 
est un polynome de degré v en z. 


Cherchons les racines de ce polynome. 


Remplaçons la première ligne par la somme de toutes les lignes. 


Chaque élément de la première ligne est alors égal à z +Y n. À 


Se 
m 


Un facteur du polynome |M | est donc z-z avec z,=- > nA==r(k 1). 
Ss=]l 
I 


Si X est le vecteur une condition nécessaire et suffisante pour que 


CPE - 


IM| = 0 est que M. X= 0 ait d'autres solutions que la solution toute nul- 
le. Cherchons les solutions non nulles et linéairement indépendantes 
de M.X = 0 ; à chacune correspondra un zéro de |M}|. Soit S.(x.) la 
somme des variables correspondant aux sèmMeS associées de la ième 
variété. Les v équations M. X = 0 s'écrivent : 


PR FD RENRS CE) 0 (=12, ...,v). 
SE 


Faisons la somme de ces équations par rapport à s: 


3Ë SE 1 
RÉPARER pren. À] S1(x:) (i) 
Proc ese 
+[(z+ AP, HAApENFEE +À, PA “0, À.]S.(x;) 
ones 


per Àp ++ n, A 1S, (x) 0 
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Posant le coefficient de S;(x;) égal à a.,, le système d'équations (i) 
shéCrit 


(ANSE SUR) ER ae SCD = OMS ENS RE 


Or la somme des v équations M. X = 0 s'écrit : 


v m 

\t Fe, = 
Ÿ x;=0 ouencore x;+ YS.(x;)=0. 
i=1 s=] 


Nous pouvons toujours supposer x; # 0, donc les S.(x;) non tous 
nuls, ce qui est vrai si et seulement si : 


HAara ti 0 (sur en) 


THEOREME 


(z-2,) (2-2) 7... (2-24) 


= r(ks1);2,,--.:21(tem)zéros distincts de | Ales 


Merle 


2/ - Cas du bloc incomplet partiellement équilibré a deux classes 
associées, 


|A | s'écrit : 
il 1 2 2 
re Ë +Ap,, “ À, P, “A,n, À Par Az Par Ai 
1 2 
ÀspaAe pe nn z + Ab A D An: 


ou 


lA=z+ (CAL AP, POS (AEAN)IAECAONNUNRE A pl 


Posons : 
Y = 2 à l 
pe PL, 
2 JL 
B d 11 h PB; 
A =Y"+2p+1 


Les racines de l'équation | A| = 0 sont : 


Z3 = 1/2{(AS A) YT CU NMNAT+X +, 
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A'estpositif, "donczf z;et/si À)> A2 0, z,< z,. 
Dans le cas où b <v, alors z,=r. 
Déterminons les exposants a;. |M|=(z-2,) (z-2,)*" (2-2 )*? avec : 
CS NS 1 
|M| peut sécrire : 
EE (2 rar oz Er REZ) (Cr) 0 CZ): 


Le coefficient de z'-1 étant nul, d'après le calcul du déterminant |M |: 
ZT EU 7 m0 L tx, =.v-1 d'où 


- {v=1) LC) y +VA+ 1) - 2n; 


Œ : 


2 VA 


Il est nécessaire que l'un des nombres &; soit un entier positif infé- 
rieur à v-1; l'autre l'est alors également, 


Remarquons que, dans le cas b < v, la relation z, = r peut s'écrire : 


pl 2 = 
HR AU ES À, )P,, (ras) P,, ]=(r- À ) (r-À,). Avec b>v-a,=a,+ 1 
Dans le cas général z, <r s'écrit : 


(A) FA ZA) Cr À) pa (r- À) p°] <(r-À) (r-À,) 


résultat dû à Nair (1943). 


Remarque. 


Le théorème dû à Connor et Clatworthy, est une généralisation 
du théorème dû à Schutzenberger. 


B - CONDITIONS D'EXISTENCE D'UN BLOC INCOMPLET PAR- 
TIELLEMENT EQUILIBRE A DEUX CLASSES ASSOCIEES 
DANS LE CAS OÙ b «v - 


Nous allons chercher tous les plans d'expérience de cette caté- 
gorie lorsque le nombre r des répétitions est fixé. 
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Comme les À, sont des nombres entiers inférieurs ou égaux àr, 
ces paramètres ne peuvent prendre qu'un nombre.fini de valeurs. 


Nous allons déterminer les paramètres de première et de deu- 
xième espèce autres que À; et r. Remarquons qu'il suffit de connaître 
be D pour connaître les autres paramètres de seconde espèce, étant 
donné le système de relations : 


np} = ERRE 
1 3e = 2 2e 
DAS tres DIS CRI ARE 
É 1 2 2 = 
= + + = 
P; ui P P,; P,; 1 Pi 


Soient À, et À, À.>À,, deux entiers positifs inférieurs à r. Les coef- 
ficients de première espèce vérifient le système de relations : 


vr = bk 
Ve Ein in er dl 
ÀAin,+A,n, = r(k-1) 


1/-Cas où À, et À, sont premiers entre eux. 


D'après l'identité de Bezout, il existe deux nombres a, et a, en- 
tiers, inférieurs en valeur absolue respectivement à À, et À,, tels que 
AA, Fa À ile 
BOSOnSRUE KES 

de a,À,ru +a,À,ru = ru 


A ni+ NA -n;=7ru 


nous tirons : (n,-a,ru)À ,+{(n,-a,ru) À,= 0 et, À, et À, étant premiers 
entreeux.: 


ou P désigne un paramètre à valeurs entières. 


a) Recherche de valeurs de b susceptibles de fournir un 
plan d'expérience. 


n,et n, sont des entiers positifs, ce qui fournit la limitation pour P: 
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a2Tu pe airu 
À; À à 


Le nombre de blocs b est donné par vr = bk, c'est-à-dire : 


HA (Rene dl}r s (a, +a,) ru + (À,-À,)+ 1 
DURE UE 


La limitation imposée à P fournit une limitation pour b et les valeurs 
de b susceptibles de fournir un bloc incomplet partiellement équilibré 
sont alors connues, b devant être pris entre deux entiers positifs. 


b) Existence de plan d'expérience de paramètre donnés r, 
HAN 2 
Soit b une valeur convenable. D'après le paragraphe précédent, 


pour qu'il existe un bloc incomplet partiellement équilibré de paramè- 
tres rb, À, À, (b<v);ilifautietiil suffit que : 


MON (EAN CNE AD) ET À )p Ar À) p llétantivérifiée avec 
P,, = t entier positif ou nul, 


_ (v-1)[-Y+NA+1]-2n, 


à 2\n 


(04 


nor 
. = nr ES de “ 
soit entier positif. De P, cut: AR Fi , nous tirons : 
n 
2 


(ONE EE EE) UE re) (ar nn) 


CREME AD ROSE] 


soit t/n, = 
(A=A,) [(r-A,)n,+(r- À,)n,] 


Désignons par N et D respectivement le numérateur et le dénomina- 
teur de la fraction du second membre, Compte tenu de a À rai ASSRLE 
nous obtenons : 


Nr air (l ra UP AC) APTE À) SAS] 
Dit) )ran(ashai)h pt X:-tx,) ru] 
Or, bk = vr nous donne : 


blutl)æ=r(ru(s +a)"F 0 (AXE X.)+ 1] 


d'ou : 
Pr(Ai-À,) = b(ur1)-r/u(ait a,)-r 
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Par suite, N = (r-XA,) [r(l+a;u(XA;-A,))- À,/r. (b(u+1)-r?u(a+a,))] 
Nr=(r-À,)(r°-X,b) (u+1) 


Dé em AS AS) (DEr OUT) 


Finalement : 


_ (= À) (r2-2,b) 
r( Ar A) (b-7) 


t/n, 


expression indépendante de n donc de k. 


Remarquons que r est toujours inférieur à b, donc que b-r est dif- 
férent de 0. 


Su) (NA-Y+1)-2n, 
2Va 


Dans le cas envisagé (b <v) 2z,=2r=À,+À,+(A,-A,)(VA-Y) 


Calcul de 


2TIAUATE A) 
À 17 Ào 
_ [ru (aita,) +P(A;=A,)]. (r- A J/(ALA,)- (a;ru-x,p) 
Samurai ein e VA Pie qu 
ŒAArSAL) SE ruiir (aa) bb (ut) TE (are) 


(b-r) (u+1) 


Gines Aie ne 


il 


Va - Y e ACQAT À,) 
tte 
AVECAY = nSn)/n. tt a tru À Pet in outattà ) Ab 
b(u+i) -ru(a;+a,)-r _ (r-;)(r2-A,b) [b(u+l)-r] 


V= a ru -— À 
: É r(À-À,) r2(A,- À.) (b-r) 


Pr SE TNA) AIR DIE RP (DER) bjr En TARN )E Ah] 
r?(A;- À.) (b-r) 

donc : 

V3- Ur br é QAR AS BEA AT RIDE Tr) END) 
r2(À,-À) (b-r) 


et, enfin : 
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(bn) A (CEE) 
u[r-r'#+br? (AÀ.) = b2 À, À;]+.(b-r) (r3-X,Àb) 


% = 


et, en tenant compte de u = k - 1, 


r?(b-r)? k 
DOTE A) (RASE (TEEN D) (FE) 


expression indépendante de a,et a,. 


2/ - Cas général : À, et À, quelconques. 

Soit d leur plus grand commun diviseur À,= 4H, d ; À,= Hd. 
D'après l'identité de Bezout, il existe a!, et a! entiers inférieurs en 
valeur absolue respectivement àh, et H, tels que H,a!+H a! = 1; de 
plus, d(u,n,+h,n,) = ru, donc d divise ru. 


Nous devons donc distinguer trois cas : 


a)AATMENTI: 
b)hu=r dus 
ChRdM= RP Avec pri Us AU, 
ue 
Nous pouvons écrire un, + bn, = ru! 
PAUL 


h, et, étant premiers entre eux, nous posons dans tous les cas : 


= Aa! Est ! 
n,5a ru b, P 


n 


tyst ! 
ra, Than P 


où P' est un paramètre à valeurs entières. Si p = P'/d, en posant 
a: = da;, nous obtenons : 


où P n'est plus un paramètre entier, mais un entier divisé par d. Nous 
obtenons donc les mêmes expressions pour t/n, et que dans le cas pré- 
cédent. 
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3/ - Conclusion. 


Le résultat remarquable est que & est toujours une fonction ho- 
mographique de k, et non comme nous pouvions le penser le quotient 
de deux trinomes du second degré en k. Les points à coordonnées en- 
tières de l'hyperbole d'équation à = ® (k) peuvent s'obtenir très faci- 
lement, Nous sommes alors en possession de toutes les valeurs de k 
susceptibles de fournir un plan d'expérience pour chaque valeur de b 
convenable, compte tenu de b > a+ 1, 


Il suffit, alors, pour conclure sur l'existence du plan d'expé- 
rience de paramètres r, À,, À,, b, k, de calculer d qui doit être un 
nombre entier et t qui doit être un entier positif ou nul. 


Nous pouvons aussi très facilement affirmer l'existence ou la 
non existence de blocs incomplets partiellement équilibrés ayant pour 
paramètres r, À,, À, (b < vjet, lorsqu'ils existent, calculer tous les 
paramètres de première et deuxième espèce. 


4/ - Cas particuliers. 
HT 
Les valeurs de b convenables sont telles que b>À,>À, 


à = - A); 
NA, D 


_ 
(l 
[æ] 


Le numérateur de est positif quand le dénominateur est négatif, 
c'est-à-dire quand À X,b. Il n'existe aucun plan d'expérience de 
paramètres r = À,, À,, b. 


THEOREME - 


I1 n'existe aucun bloc incomplet partiellement équilibré ayant 
pour paramètres b, À, =r, À,<À,(b<v)si À?<A. b, 


Nous allons dresser une liste des valeurs des paramètres À.,b et 
r<10 pour lesquelles il n'existe aucun plan d'expériences. 
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2) <eD. 


Dans ces conditions t est toujours négatif, Nous pouvons donc 
conclure que quel que soit À, il n'existe aucun bloc incomplet partiel- 
lement équilibré de paramètres b, r, À, (b<v) si r?<A, b. 


Le tableau ci-dessus est donc vrai quel que soit À,> À,. 


3°) Remarque : 
Dans lecas r=À,, nous retrouvons un théorème connu, à savoir 


que tout bloc incomplet partiellement équilibré de cette espèce s'obtient 
à partir d'un bloc équilibré en y remplaçant chaque variété par un 


groupe de n variétés, 


En effet : 
&œ= —= —— d'après la définition même 
Posons : 
PP NID NV EE) EC) 
r(k= 1) =rn; ce x(n il) 


Ée (En) (rrax) 
1 


1 
D 


Posons : ds 


1 + 


2 
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rm 


Uor+ (=D) 


Nous déduisons de l'existence du bloc incomplet partiellement équilibré 
l'existence du bloc incomplet équilibré de paramètres : 


Réciproquement, supposons l'existence d'un bloc incomplet équilibré 
de paramètres : 


v* r* = b* k*—> pb = rm 


À (m-1 = ee 
_ = _ ——+ EE ——————…——————————— 

a unre r+A\(m-1) 

nous en déduisons donc en remplaçant chaque variété par un groupe de 

n variétés, l'existence du bloc incomplet partiellement équilibré de 

parametres: v=nmm bb, rar =k= pr, \2r, As À On TRE 


n, = n(m-1), p,, = 0. 


CHAPITRE Il 


BLOCS INCOMPLETS 
PARTIELLEMENT ÉQUILIBRÉ A 4 RÉPÉTITIONS 


Nous allons énumérer tous les plans d'expérience de cette caté- 
gorie. Pour cela nous résoudrons le système d'équations diophantien- 
nes, imposé aux paramètres de première et deuxième espèce. 


vr = bk 
1 VAE EEEr Cri 


Any, n, =r(k-1) 
Avec mrimb.an.>"0,entiers 14 > À: > 10entiers, pé = n ed. 
n 
F2 


avec pe > 0 entiers. 


Le nombre & déjà utilisé dans le chapitre précédent étant entier 
positif, 


Nous savons de plus que tous les paramètres doivent vérifier 
l'inégalité : 


@) (r-A) (r-A,)> CAS A) Lr- À) p£ - (r-A.) p£,] 


Pour résoudre ce problème nous allons considérer les diffé- 
rentes valeurs que peut prendre le paramètre k : il est obligatoirement 
supérieur à 1 (chaque bloc ayant au moins 2 parcelles). 


KÉ=N2 KES KE alors pv: 
k >4, b <v nous pourrons appliquer les résultats du chapitre 


précédent. 


Les nombres À; étant obligatoirement inférieurs à k, nous exa- 
minerons toutes les valeurs que peuvent prendre les À pour chaque 
valeurs de Kk. 
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R.C. Bose, W.H. Clatworthy et S.S. Shirkande ont classé les 
blocs incomplets partiellement équilibrés connus en 5 types suivant les 
schémas d'association (c'est-à-dire suivant la façon dont sont formées 
les classes d'association), nous garderons leurs notations. 


! 

1/ - Groupe divisible (GD). 

Les variétés v = nm sont rangées en rectangle une variété v; 
est première associée des n-1 variétés de la même rangée, et deuxième 
associée des n(m-1) autres. 


2/ - Simples (SL). 
Chaque fois que À,est nul. Le plan d'expérience sert lui-même 
de schéma d'association, 


3/ - Triangulaire (T). 
_ n(n-1) : 
Le nombre v est de la forme v = rer les variétés sont ran- 


gées en rectangle de telle sorte que : 


a) la diagonale principale soit vide. 


b) les variétés soient rangées au-dessus de la diagonale 
principale. 


c) le dessous de la diagonale principale est symétrique du 
du dessus, 


Pour une variété v: les premières associées sont sur la même 
rangée (ou la même colonne). 


NS NE 


n, = ; (n-2) (n-3) 


4/ - Carré latin (LS). 
v=n?, nous rangeons les variétés en carré la variété v, est pre- 
mière associée de celles qui se trouvent, soit dans la même rangée, 
soit dans la même colonne. 
mn 2(n-1) 


n, = (n-1)? 


5/ - Carrés latins orthogonaux. 


Nous superposons au carré des n? variétés, un carré latin ortho- 
gonal formé par les n lettres A, B, ... N. 
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Chaque variété sera lère associée des variétés de la même ran- 
gée, de la même colonne, et des variétés correspondant aux mêmes 
lettres. 


n = 3(n-1) 


(n-1) (n-2) 


[l 


n, 


6/ -Cycliques (C). 


Les premières associées de la variété i, s'obtiennent en ajoutant 
i-1 (modulo v) aux variétés premières associées de la variété 1. 


En faisant la classification des plans d'expérience envisagés 
(r=4), il nous est apparu qu'il existait d'autres types de schéma d'as- 
sociation que ceux cités. 


Nous les considérons comme une généralisation des types carrés 
latins : 
1) v=nm, v; sera première associée avec les variétés de 
FX Ps Fe à ZE 
la même rangée et de la même colonne 


fe = Do un 10 er PA 

n, = (n-1) (m-1) 
2°) v = nm, nous superposons un rectangle formé par les 
permutations circulaires des lettres À, ..., N ou À, ..., M et v, sera 


première associée des variétés de la même rangée, de la même colon- 
ne, et des variétés correspondant à la même lettre. 


nil) n-lnpout les lettres À,2..5eN 
n (nl) Mm-I10 pour les lettres A, ..,1M 


I - BLOCS INCOMPLETS PARTIELLEMENT EQUILIBRES, ou r=4, 
et k=2- 


4v =2b 
2 EE 
NA EM tt nl 


La résolution de l'équation diophantienne À,n, +À,n, = k nous 
donnera les valeurs de n, et n.. 
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Nous pouvons ainsi calculer v et b. 


Puis posant pa , pÆ =n, -t devant être positif ou nul, nous en 
déduirons les valeurs Fe { convenables, elle devront de plus rendre « 
entier positif et satisfaire à A. 


Ici 2n,+n, = 4 a pour seule solution n, = 1, n, = 2 donc : 


b=2v=8 
NN ee t doit être pair 
LA ne 2 pair 
(A) s'écrit 6 >3t - 2p°, d'où t£< 2 
SPC AU 
ï P,, P, l 
SE \ & 
p P;, < Di, 
pour t = 0, 


L=VY" F2B+ 1-4 


> L( (v-1)(VA-y+1)-2n,]=1 


la valeur t=0 convient. Pour t=2, Y= -1, B=3, A = 8, « est irration- 
nel, cette valeur de t ne convient pas. 


THEOREME I, - 


La classe des blocs incomplets partiellement équilibrés à deux 


classes associées, oùk=2,r À,= 2, À2= 1, contient le seul 
plan d'expérience TM DELLE Pnm2 P,, = 0, 


EEE TE 


| 
2n, = 4, posons n,=X | 
| 


vs 34% bD'=26+%X) 


p? =2- doit être entier positif ou nul, 


donc ME RIEZ 


soit t=0 
X=t 
X = 2t 
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l'inégalité (A) s'écrit 4>4t - 2p2 


pour t=0 (à) est vérifiée, «= = est entier positif si 


hu entier, t = O0 convient 
pour t = 0. (A) donne PA diouree 

a = 2, t = X = 1, est une valeur convenable 
pour X = 2t (A) s'écrit ATLAS, MSOIATE= PM SOIT RO 


t=0 a déjà été étudié. 


t = 1, «à = 2 valeur convenable, 


pour X=2 @) SUÉCRICALES 


wo 


,; soit t = 0, 1, cas déjà étudiés. 


THEOREME l, - 


La classe des blocs incomplets partiellement équilibrés, où 
r=4, k=2, À,=.2,, À,=0, contient : 


Pile Diandierpériencenm = 2 ne 1 vend" D=T6 TE "TT, 


HP VSD, DECITRE 


l 
o 
5 


2) le plan d'expérience n, = 


8) la classe nn; = 2, n,= 34, v=3(u+1), 


n, = 4, posons comme précédemment n, = X 


2 =4 00 donc soit ù 
Pre KES X 


[I 
PI 

Ë 
N 
+ 
& 
+ 


3 
() s'écrit tU+SYE 6 
t = O donne : 


a+, donc = SU; lrentier 
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t = 0, impossible car X serait nul. 


t = 1, ä = 3 valeur convenable, 


DOME 
! 14 
= (0 4 = — 
ta=02; 3 
t=3, à = 6 valeur convenable. 
t = 0, a déjà été étudié. 
t = 1, À = 8 est irrationnel. 
X=2t, es t = 2, à = 4 valeur convenable. 
t = 3, = 12, & est irrationnel. 
ROUE ICE M PE ITNONT L 
t = 0, a déjà été étudié, 
t= 1, = 13, & est irrationnel 
t = 2, = 12, 1 11 11 
NEMATLS “ 
t = 3 = Isis nt "1 11 
t = 4, = 14, valeur convenable, 
t=5, = 21, à est irrationnel. 


THEOREME LL, - 


La classe des blocs incomplets partiellement équilibrés où 
m2, FANS l/,\ = 0;1Contient. 


1/ Les quatre plans d'expérience 


2/ La classe de plans d'expérience : 


H entier n,=4, n,=5, v=5(u+1), b=10(u+1), por il 
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II - PLANS D'EXPERIENCE ou r = 4, et k = 3. 
LINE 2, À, 2 3n,+2n,=8 


Les seuls entiers positifs vérifiant cette égalité sont n, = 2, 


n, = 1, d'où v = 4 et 4v = 3b, b = _. c'est impossible. 
PURE 9; Noel 3n,+n,= 8 


deux solutions : 


Mere: * PRAUS s 
nn. 2, V= 5, b'=—— "impossible 


3 
. = Le 207 ; 
he ils = Gi Ne, (0 = impossible 
3/ NAS TN A0 INT END, 


cette équation n'a aucune solution entière. 


THEOREME Il, - 


Iln'existe aucun bloc incomplet partiellement équilibré ou k = 3, 
DÉEMN2S, À 150: 


DS 0 VIER0 Den IeSstipas entier 
Hot cos die robe ou 11 nt 
HS IL 2 VAE 0 b=8, (A)s'écrit t<2 


t = 0, «= ; impossible 


t =. 1. pa est-fractionnaire 
L2 


t = 2, aäœ= 2 valeur convenable. 
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THEOREME IL - : 


La classe des blocs incomplets partiellement équilibrés où 
r=4,k=3, À,=2, À,= 1 ne contient que le plan d'expérience 


b = + (5+ X), une condition nécessaire d'existence est donc que 
+ — 


v=5+X=3u, L entier positif. 
(A) s'écrit 2>2t- (4 - t)Ravecs 
4 3h -5 
soit t = 0 


SOL 3 UT 51 2, dite 


Or t >0 implique u>2 


la seule valeur convenable est p=2,t=1, a=3 
c)3u=5=2t, (A)s'écrititez 


t = 1, irrationnel 


t=2, «a =4, H = 3 valeur convenable 
GRGUMONSREN (A) s'écrit tea donc t = 1, 2 


t = 1, 2 œest irrationnel. 
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THEOREME Il, - 

La classe des blocs incomplets partiellement équilibrés où 
PS4 Ki=:3; À,=10 contient 

1/ la classe n,=4, n,=5(3®%- 1), v =15%9, b= 209, t 
® entier positif. 


2/ le plan d'expérience n,=4, n,=1, v=6, b= 8, t = 1 note 
SR, dans la classification de Bose. 


3/ lepland'expériencen, =n,=4, v 12, t = 2 note LS, 
dans la classification de Bose, 


et pour plan : 
C2) 2006) 
(4, 5, 6) (4, 5, 3) 
(1, 5, 6) 
(4, 2, 3) 


Le plan d'expérience LS, a pour schéma d'association : 


Su 
© 1 N 
© OO © 


et pour plan : 
(1, 2, 3) (1, 4, 7) 
(PS5 R0)I(272578) 
(20/29) (576709) 
répétés deux fois. 


GRAS EL = 10 n, = 8 posons n,= X 


— b entier exige 9+X=3k 
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i iti 2? = 8-———t implique 
U entier positif, et P!, 8 3(h- 3) pliq 


soit t=0 
SO AIN -5)Et, 21, 40 61 


SIT 


a)t=0 (A) s'écrit 12> - 3 x 8 toujours vérifié a=—- 1 
qui exige 1 = 39 , ® entier > 1 
d'ouan,.r 8, 1h,1= 991) =v2,90, De 12 st 0 
bj3Qu3)k=" t (à) s'écrit u<4, or t >0 exige U > 3 
u = 3, X = 0 impossible 


uU =4, t = 3, & = 9 valeur convenable. 


c) 3(U-3)=2t (A) s'écrit t < 6, t devant être divisible par 
SU MtU=NS QUE At I6 COR IrAaCtiOonnarre 


1 


Oo 


d) 3(H-3) = 4t (A)s'écrit t <=, t multiple de 3 


| 


t=3, = 7, ÀA=28, à est irrationnel 


t=6, d = 11, À = 25, & est fractionnaire 
e) 3(u- 3) = 8t (A) s'écrit t &  t multiple de 3 
D 4 : 


t= 3, nu ="11,A= 37 "a est irrationnel 
t=6, = 19, = 28, 1 11 


THEOREME Il,-- 


La classe des blocs incomplets partiellement équilibrés où 
TASI4 KES ÀL= I, "ASE 0rcomprendi: 


1) le plan d'expérience note SR,, dans la classification de Bose: 


= 8, M re 12, Del it 


2) le plan note T;:n, = CNET NN AEMEE de PA), ie 


3) la classe de plans d'expérience : 
NERO, 0 RONDE) VA RO ID DEN ET 


ou & est entier > 1. 
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SR,, a pour schéma d'association : 


et pour plan : 


(1, 2, 
(1203, 
( 9, 10, 
(E1,64; 


(51,268, 
( 7, 6, 


(AAC T2, 


6) 
4) 
7) 


1407, 10 
CRE EI 

RG 0 T2 

MCE 2 y DL (LA 5: 
7) (5, 7, 3) (8, 3, 
HACEURNCRERE 
60 2200)"(27.127 10) 
9) 

2) 

3) 

5) 


T;; a pour schéma d'association : 


et pour plan : 


2 


© À J © © OO ON 


© O1 © © 1 À © © — 
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III - BLOCS INCOMPLETS PARTIELLEMENT EQUILIBRES ou k=4 - 


1/ °À;="4 (a) s'écrit t £0 la seule solution possible est donc 


LEO 


a) À,= 3 l'équation 4n, + 3n, = 12 n'a aucune solution en- 
tière positive. 


b) À,=2 2n;+n,=6 a pour solutions 
soit -n, = 2,-n, = 2, avec a = Squine convient pas, 
soit n, = 1, n,= 4, p}, = 1, «= 2 valeur conve- 
nable 
CAS =T Ian tin #12 
SON MN, 20, Ad ONE 
soit n,= 2, n,=4, 4=9, a fractionnaire 
d)À,=0 n, =3 posons n,=X 


PRG, AE a= À, d'ouX=4p, uentier > O. 


THEOREME III, - 


I1 n'existe aucun bloc incomplet partiellement équilibrés ou 
CR AR EN NN Sd Un 


Il n'existe qu'un seul plan d'expérience où r=4, k 
N, =2, il est notés 


ns nn 


Iln'existe qu'un seul plan d'expérience où r=4, k= 
= d/ilestnotess; 


n,=1, n,= 8, v = 10, b= 10, t=0 


Il existe une classe de plans d'expérience our =4, k=4, 
À = 4, NE 0, 


1,n,=4u, v =2(2u+1)=b, t = 0 ou u est entier positif 


S, à pour schéma d'association : 
IN ES 
4 5 6 
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et pour plan : 


ST 25 25 detre 
2 OSIG D ZONE 
SONIA CEST 
S:, à pour schéma d'association : 
INDE AS 05 
GAMME 
et pour plan : 
D 10 1 lg & À Z 
on Gr Loan CSS 
Po DU EE ton 200 9774 
SCT 9 CS RAIIOSS 
4 © & 0 CC MT 
AIX =, DA T2 3n, + 2n, = 12 a pour seule solution 


12 


(à) s'écrit 2>2t-p° @+i)t-2 


p, entier positif exige t = 3 pavec = 0 ou ®=1 


(a) s'écrit alors ®< J. donc seul ®@ = 0 convient A=9, «= 1. 


THEOREME III, - 


Ilexiste un seul bloc incomplet partiellement équilibré où r =4, 
ne ne | me 


ilestnoté R, dans la classification de Bose:n, = 2, n,=3, b=b=6, 
P 


12 


C'estun''Group Divisible'', son schéma d'association est : 


ll 6 © 
2 4 6 
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son plan : 
102,4 6 
SJ CE 7 
SG 2 € 
IDE SOS 
LS 
CO PES 
DRASS A El 3n, +n, = 12 cette équation admet 3 cou- 


ples de solutions : 


anne n59; pe : les seules valeurs convenables 


de t sont donc 0 et 9 
@ ; ; 45 ; M Û _ à ; 
s'écrit t<ss ce qui exclut la valeur t=9;1t=0, «= 2 et n'est pas 
entier. 


b)n=12;urnu=265 Ip" = 2-1 d'oùt=0, 3, 6 


12 


(a) s'écrit t< . , ce qui exclut les valeurs 3 et 6 t=0, A=9, a=2 


chnr=3;n DD eine d ON SEIeRO: lu 0 


1 


(a) s'écrit t < è ce qui exclut les valeurs 2 et 3 


DID EI e n'est pas entier 


t= 1, A= 8, « est irrationnel 


THEOREME III, - 


Ilexiste un seul bloc incomplet partiellement équilibré où r 
REPAS ES EX 


2 


ilestnoté R, dans la classification de Bose :n, = 2, n, = 6, v=b 
P 


12 


C'est un ‘'Group Divisible'', son schéma d'association est : 
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D er 
2 
36.9 
Son plan : 
PAT? 
ASC ES 
JOIE 
A Ts ln 
Ge) 276 
CMOS 
LA 
22/0509 
ei 6) Gp Al 
n'est posons nn x 
e -4-% ce qui exiget = OouX=t, 2t, 4t 


: at 
= ! - — ifi = 
a)it =00 @)s écrit 1 > 4t - 4 + x donc est vérifiée A= 25, 
a = d'où X = 51 avec L entier > 0. 
b}xX=t (A) donne t <1 d'où t=X=0 ce qui est impossible, 
ChECEI2Tt (A) s'écrit E & + doute A0 


t = 1, À = 8, & est irrationnel 


t = 0 a été étudié. 


5 
= ' : NES = 
d) X at (A) s'écrit t <4 d'où t 1e «0 


É 


1, ÀA= 13, & est irrationnel 


0 a été étudié, 


t 


THEOREME Ill, - 


Les blocs incomplets partiellement équilibrés oùr = 4, k=4, 
A1= 3, À,= 0 sont formés par la classe des plans d'expérience : 


n,=4,n,= 5h, v=b=5(u+1), Pere 0 b entier > O0 
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l'équation 2n, + n, = 12 admet 5 couples 
de solutions. 


t 
00) CO OO 


a), € 1 & 


@) s'écritt< . ce qui exclut la valeur 10 t=0, A=4, à =5. 
b = 2 = 8, p,=2-+ d'out=0, 4, 8 
) n É. » M2 ; pre 4 u , , 


(a) SIÉCRLURTE F ce qui exclut les valeurs 4 et 8 


t=0, A=9, a-= = donc ne convient pas 


t » 
Css ne s16; bare be 0, 2, 4,26 
@) SHÉCTILE TES LRO NE : ne convient pas 


LMD NAS O0 C5 
= = 2 = _. = 
din te Up FAT DUFAO RE 22, 03,14 


1 


(a) s'écrit t <É ce qui exclut t = 3, 4 


E2=20; 125, AU . ne convient pas 


t = 1, A= 13, à est irrationnel 


t=2, A == = A 
CA 215,0 22 D = 002 


(A) s'écrit t<2 d'où t=0, 2 


t=10 05-1636, die 
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THEOREME IL, - 


Les seuls blocs incomplets partiellement équilibrés ou r=k=4, 


ÀA1= 2, À,= 1 sont les plans d'expérience de la classification de 
Bose : 


R:; a pour schéma d'association : 


Es Ki nu TS 
DK ©) HO) SU 2 


Ev Dour Jen nt, be Hbc MANETTES RE nil 


T,, son schéma d'association est : 


A ET 
TS ENG 

2. 5 & GS 9 

36 LO 

ART OO 

son plan : 

7 do € 7 
HO Mn 2 5 
TRS Ce? 
Gen 2 RP NX! 
1 9 10 8 
Hd 3 210 
ST NL ET 
TOR ENT 9 
JAAGr 1 3 
CR 5) 
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LS, son schéma d'association est : 


1726) 

4 5 6 

US, © 

son plan : 

ne LC) 7 
CRETE 
SOTLNEA 9 
C2 SET 
2 OS ATERS 
AN OBS) 
5% 106 
ue Gr & 
315, 8 1 

GENE 27 À =) 0 n,=6 posons n,=X, p?,=6(1- +) 


d'oùt 0m ou Et 2 SE Gt 


a)t = 0 (A Jest satisfaite A=(7) a= À donc X = 7H Hen- 
tier > 0 

b) X = t(A )s'écrit X<1 donc X=1, A=4, a=4 

c)X = 2t (A )s' écrit t< 2: te L 2: 


= 2, A =S192 ain 1 


| =1, A=13 «est irrationnel 
d) X = 3 (A )s'écrit ts 3:t=1, 2, 3, 


À = 20 «à est irrationnel 


15 
2: A=17 œ On" nt 
3, 


U 
t 
t 


A=16, a =6 


e) X = 6t(A s'écrit t «hits 12,305 
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t = 1, À = 29 à est irrationnel 
END AE240 at "1 
= AE IE n 


THEOREME III, - 


Les seuls Blocs incomplets partiellement équilibrés où r=4, 
À,= 0, sont : 


1) La classe n,=6, n,=7H, v=b=7(u+1), PRea0 


L 
12 


3) = sv 10 D 
al 
4) 9, v 16be 


2) Le plan d'expérience SR,!n, =,6, n,=.1, v«= b=.8,. p 


SR, a pour schéma d'association : 


25 nt 


et pour plan : 


© BB J © D NN Où mn 
En mm æ œ@  J D NN 
D OO OO mm æ OO 
D NI NN D OO mr œo À 


T, a pour schéma d'association : 


x 2 T7 
TES 50 
LS ER OURS 
A Gr aiG Se Ait 
1 & €) de 
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et pour plan : 


101203 mr 4 DETTE 

5 G 1 1 GES IT 

SOS OR 52 

1e) Er € 8 STORONC 

ANTON OEIL 0 T4 10280 
TIRER 10 n, = 12, posons n, = X 


t (5 Et, 
DE _— ‘Où 
pi, = 13(1-3) d'où ee =t, 2t, 3t, 4t, 6t, 12t 


a) 0 (A) est vérifiée a =, donc X=13hH, entier > 0 
b) X =t (A)s'écrit X < 3:3: 
Mist ls Dean ar 


X =t=2, A=9, &œ = 10 
RITES, ANG RER 


c) X = 2t (A) s'écrit t£e= d'où t = 1, 2, 3. 4, 5, 6,07 


t = 1, À = 40, à est irrationnel ; il en est de même pour 


teN2, 3000, 1 


Se n'est pas entier 
10 P : 


LIGUES? 5 RTE 


LE A 


d) x =)3t (A)stécrit.t <9, d'ont= 1. 23/4/5078 


TS 2 SA 5 06 LT est irrationnel 
dE 6h Me MOI 
t = 9, À = 36, x = 24 

ex rat (à) s'écrit te © d'où t=1, 12. 09 
te 2 3 4,05; 16,170 est irrationnel 
t=L8 A = 36, x" 20 
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f) X = 6t (A) s'écrit t< ©: d'ou te 1 0 10 


Fe li2 03 & SAT 2, OLI0, 
t=6, A=49, «= 12 


est irrationnel 


g) X = 12t (A) s'écrit << d'ou Pets sil 
est toujours irrationnel. 
THEOREME III, - 


Les seuls blocs incomplets partiellement équilibrés, 
EMI =X, À) =01, .À= 0 sont: 


1) la classe des plans d'expérience 
Hrentier. >0 


120 13H, = pbs 18 (#1) p} 0 


2) les 4 plans d'expérience de la classification de Bose : 


3) les 4 nouveaux plans d'expérience : 


v = b = 25, 


45, 


R,, à pour schéma d'association (groupe divisible) : 


et pour plan : 


253 


où 
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2 3 5 8 LORS DRICR? 
3 4 6 9 US AN ASS 
4 5 RRRLO 14 SMIIONR6E 
5) 6 Ali 13 14 SES 
6 7 22 8 (Jon CS Cr 
Ÿ 1 nil 12 13 8 4 
1 2 4 14 

Ge some 1 


R,, ‘'Group divisible" a pour schéma d'association : 


et pour plan : 
on Es Dern 
module 15 


AVEC DIE Id 


SR, ‘'Group divisible"" a pour schéma d'association : 


Fo © de 

2 6 10 14 

Set ile 16 

AO 2 LC 
et pour plan : 

2 3 1 4 LU ATEN € UN 15 
LS 10 8 9 6 9 12 3 
lg Ait A KE 4 13 T0 
14 12 5 ti 8 5 2 il 

tu Se 6 1 10 LATIN? 

3 16 10 5 9 7 16 2 


il 4 9 14 13 14 3 8 
12 CLS TS 5 6 4 15 


CONSTRUCTIONS DE BLOCS INCOMPLETS — ANALYSE DE VARIANCE 255 


Si 2 

5 6 « 1 3 Ole g OSSI 26 
15 ni NO 410? 5 OS 62710 
23 24 25 CN SOMS CR? 2ONMOGN25 2 
ill. 66r ir 07 14 C2 OT 31. 23 Si 
2 8 OL 18 SP 11 30.78 4 23 
ZT Se 6 19 14 21 4 CCM VPN 
CRU 5 25 SOS CS O0 LEE TONI? 
26 19 6.36 12 Go 1 Gil SOS JOINT 
20 CSSS Lo Où 29 JE at 22 25 8 
CLS RS) ec (0 7 816-2840 2 02/1 OS 
40 “36 - 32 3 32 ASIN 3T 6e. Zn I ZX 
32 4 39 37 Sem nnilile, 2 DORE ESSS 
2820."33 15 34 31 e) © 17 OMUSO U3A 


107 34 18 29 


IV - BLOCS INCOMPLETS PARTIELLEMENT EQUILIBRES ou[k >4.] 


Danc ce cas k > r, donc v > b, nous pouvons appliquer les ré- 
sultats du chapitre précédent. 


Nous avons vu que dans ce cas 
E & (0) @t 


THEOREME IV - 


Les seuls Blocs Incomplets Partiellement Equilibrés où k > 4, 


À,=4, À, #0, sont obtenus à partir des Blocs incomplets équilibrés 
ayant pour paramètres: v”=m, b', r°=4, k* =p, À = 1, 2, 3, en 
y remplaçant chaque variété par un groupe de n vanetes, 


b- 


ZAR TE NES & =2=#, donc b=4(1+a),n,=k-1, n,=ak 
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THEÉOREME IV, - 


Les seuls Blocs Incomplets Partiellement Equilibrés où k>4, 
À,=4, À, = 0 forment la classe de plans d'expérience : 


=ak, v=k(l1+a), b=4(1+a) p,, = 0 


où à est entier positif. 


JJR Ni Se +2 Nous reprenons les notations du chapitre 


précédent. 
3n, + 2n, = 4u 
etr: el LES dobea.f tata, di 
n, = 4u.- 2.0 
> MER AUSTIN n,> d. donne, ÿ <2u 1 
n > lRdonne p « 
= AU PE) : b(u+ 1) 
Lonrere Ge TT 
Cor: . 16 u 
la limitation sur © fournit: 3 < DE< Ti 
6 be 
L 
a)b=6 ta) avec & AuUb- 1 +5 doncu < 9 
4+7 à 
Re mr doit être entier 0, et u> 3 
War fl u + 7 doit diviser 3.2 


U +7 3/97 NE 9) DATE Det etes 


u +7 = 2" u=9", Ps ns nn, 644t=3 
3°, 2/(u+7 7u + 3 
b\ibi=27 a = 2-2"(ur7) [04 = U 
) Bu + 33 avec <OMeCtRE 2 
a a 
Vs ALES = 14+@, entier 


5  5(5u + 33) 


CR iL _5u+9 
de plus ont avec.n,= 2 
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Il faut donc que 5u + 9 = 24,X a < 6 entraîne 5u + 33 < 3.2* 
SRE De do a 
it 
14.2 


toutes ces valeurs rendent ® entier, mais &> 6 donc ne conviennent 
pas. 


THEOREME ee 


Les seuls Blocs Incomplets Partiellement Equilibrés de la 
classe r = 4, k >4, À À =12 "sont : 


LION STI (A) s'écrit 3 = 6t - 2p£, 


Cette équation n'a aucune solution entière et > O. 
THEOREME IV, = 


I1 n'existe aucun Bloc Incomplet Partiellement équilibré 
RS 4NAE= 


HA: =S, À,= 0 (à) s'écrit ADI prob 


Cette équation n'a aucune solution entière et positive. 


THEOREME IV, - 


Il n'existe aucun Bloc Incomplet Partiellement équilibré où 
REA RATES; 


> 1, Mn: donnent £i,;:,.20, <P< Auy-21 
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b(u +1 


AVEC à 


finalement 8<b < 15 


a)b=8 az + (u+1), n,=2, t=1 


a<b- 1-7 


I1 faut et suffit pour assurer l'existence du plan d'expérience que 
k=u+1= 3h 
avec 4u<7, c'est-à-dire H<1 


ON AONNENREESETT 


Il n'y a donc aucun plan d'expérience de cette catégorie avec b=8. 


E RENNES DNS 
DES, ss 8 rérmaeus) 0e 


a<fexige -prentier >.0 


a>0 5 S = < 8 


t 7 
de plus Ce 


10” 2n, = u +5 doit être multiple de 20. 
2 


La seule valeur remplissant ces conditions est : 


Tu. +,115= 22007 5Fouct 45-1760 
Dis; 1K 275617212606, n 1920 
. _ 24(u+1) 
c) b = 10 RE £'I,rexige u <5 
Le 
= Ste) u + 1 doit être pair 

,-M2- u+3 
2 2 


La seule valeur de u convenable est, en fin de calcul u = 5. 


72, 8(u+1 2 Te MOT 
Dib= ii gl 7 ren 
) 15 usa UE Mao 
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Cette condition s'écrit finalement u< 3 il n'y a donc aucune so- 
lution, 


2e L 003727 


Go © 12 CORRE E ÉROCE” 
u+10 u+10 


CHASAIOURUSS 


Il n'y a aucune solution. 


3, 2}(u+1) ARS NS GE) 
f\#D=2813 ER PTE, QE Ta PES Re is 
) But 57 drone 
Cette dernière condition s'écrit u $3 il n'y a donc aucune solu- 
tion. 
200 (u+1) 2 200. 14 
D= 14 (RE E ER EU _ ————— 
8) Suds MST 0 


Condition qui s'écrit u <7 ; or toutes les valeurs de u qui ren- 
dent & entier sont supérieures à 7. 


27, 11° (u+1 2 180 
b = 1 So LI AUET) = À ———_— 
D) RTE TEEN RC ETE RUE 


14u + 194 
= + —————————— 
138 yorur187 


or 194 n'est pas divisible par 7. 


I1 n'y a aucune solution. 


THEOREME IV, - 


Les seuls blocs incomplets partiellement équilibrés de la clas- 
se r = 4, k> 4, À,=2, AE 1 sont : 


IDD TAnET D = 10, KE 6, VOLS, td, DE, 6.0. 56 


2) le plan de paramètres : b = 9, k = 756; vw = 1701, t=266, 
HS 20: 


V - CONCLUSION - CAS D'IMPOSSIBILITE - SCHEMA DES NOU- 
VEAUX PLANS D'EXPERIENCE TROUVES PRECEDEMMENT - 


Chaque fois que les plans d'expérience obtenus, étaient déjà con- 
nus, nous avons donné la référence de la classification de Bose, ainsi 
que le schéma d'association et le tableau du plan. 
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Nous avons démontré, que dans 6 cas il n'existait aucun bloc in- 


complet partiellement équilibré de la catégorie envisagée ceci pour les 


valeurs des paramètres : 


Fed ke VON ms A2 
Led ok 9 0 AVE IONIMENS 21 
De KES AS CR EO 
D = 4 MES 4 MONET, ie 
r=4, DA 0: À,= 1 
pe Ci k > 4 A9; À,=0 


cé 1 


En recherchant systématiquement tous les blocs incomplets par- 
tiellement équilibrésoù r=4, nous avons démontré l'existence de nou- 


veaux plans d'expérience ne figurant pas dans la classification de Bose. 


Nous allons donner ici la liste de ces plans, leurs schémas 
d'association, et leurs tableaux. 
LUKE ONE 2, LXsil; 0 via4 Dr 80 #1 in =2,p 0 


L 2 


C'est un ‘'group divisible'', son schéma d'association est : 


son tableau : 


2 )Rk=E27 ATEA2 À5 = OV DO in 202 0 ND I 


2 12 


C'est un ‘group divisible'', ayant même schéma d'association 
que le premier. 


Son tableau est : 
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C'est un plan simple : 


PRO ALMA O0 Ve Gb #12, 4; n,= lp 


12 


C'est un plan triangulaire. Son schéma d'association est : 


Ka LS 2 NRA 
Lx 325) 
2 XU6 
ARS OCR 


son tableau : 
1e Ca Le 0 PET MS NS Nr er 
D LG CS 5 GG 


DARERR A Et À, D V8) beG n=4;n 3,p =. 


C'est un ‘group divisible"', Son schéma d'association : 


122 4 
SRG 7 ES 
son tableau : 
4 ai Ge ES A TT ER Ce 7 
2030 41504 4 7 TES 8 8 


Dis, À1= 1 N=10, Viey21; bi=;42; n;=#4,1n7.4106, DEEE 


C'est un plan simple : 
1 1 il TC MC EC MOT RETRO TRS 
2 3 4 DRAC SOS TPS MON IORS MOTO MIO 


SO At MES AN VA ÉME RINGEN CT t TEE TS: 
OR LS AE CS AT OO 202 O2 019 


18.19" 18.20 
20) AN AAA 
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TAKE, NE, À SE Levi, bn MER EME ER 


C'est un ‘'group divisible"", Son schéma d'association : 


12 
AN SG 
son plan : 
ES LS EE FN one er 
DELA SNS EE LR CS 
SO MOD COS NC 


BAR = ATMN = 2, À, 0 vie. 16, b = 16, n°6, 1, 19, pe 
ER A TS A hr om mp 2 


C'est un carré latin. Son schéma d'association est : 


1 2 3 4 
5 OS NET 
6 Ce mr TE 
CROIS AG 


son tableau : 


5 CN A1 0 9 8 
6 OO ST 2 
TOMOMISMIC 1601512 


répété deux fois 


BB © NN mr 


Oki 4 À = 1 (À 0 veei25, be. niela, mn =12 P,, =6. 
M ARS LT A NT me de le 
C'est un plan d'expérience type carré latin. 


Son schéma d'association est : 


H 
[En 
CNONMEMEUNT 
N 
QNONMEMEN CT 
= 
CO 
ONE A 
+ 
FN 
EN OURQORCT 
FH 
Oo 
PAU ONCE 
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son tableau : 


il 1 1 1 
812 17-22 
9 14 18 23 


LOMME 70220010 


TOR EL MEX 


d 


5 
7 
8 
9 


23 
25 


14 
15 


16 
17 


20 


4 
6 
7 
8 


10 
19 
21 
24 


13 
15 


16 


20 


C'est un plan triangulaire. Son schéma d'association est : 


1 
X 
3 
5) 
8 


12 


2 
3 


4 
5 
6 
x 


10 
14 
19 
25 


7 
8 
9 


11 
12 
13 


16 
17 
18 


22 
23 
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22 
25 


er El a NP es PP ES ras a es 


CHAPITRE IV 


ANALYSE DE LA VARIANCE D'UN BLOC INCOMPLET 
PARTIELLEMENT ÉQUILIBRÉS À m CLASSES ASSOCIÉES 


L'analyse de la variance de ‘Bloc Incomplet Equilibré'" et celle 
du ''Bloc Incomplet Partiellement Equilibré" présentent de nombreuses 
analogies provenant des analogies de leurs définitions. Nous allons les 
mettre en évidence en établissant l'analyse de la variance du Bloc In- 
complet Partiellement Equilibré dans le cas le plus général. 


Dans les deux cas, nous considèrerons le plan d'expérience com- 
me une double classification, suivant les variétés et suivant les blocs; 
ce qui nous permet d'obtenir deux sortes d'informations provenant des 
expériences : les informations ‘'intra-bloc'' fournies par la comparai- 
son des variétés à l'intérieur des blocs, et les informations ‘'inter- 
bloc" fournies par la comparaison des variétés entre les blocs. 


Nous utiliserons pour représenter le schéma du plan d'expé- 
rience, la matrice d'incidence, déjà rencontrée dans le chapitre II. 


N = (n;;) où n;;, = 1 si la ième variété apparait dans le jème bloc. 


n;; = O0 dans le cas contraire. 


Dans les cas de blocs incomplets équilibrés et partiellement 


équilibrés, les n;; sont liés par les relations : 


Ÿn;;=k nombre de parcelles d'un bloc. 
L 
Dn:;=r nt ! répétitions d'une variété. 


J 


> ni; zr 1! LL LA] [A 11 
J 


Sin A, L "" blocs contenant à la fois la variété i et la 
À variété s. 


dans le cas du bloc incomplet équilibré À,.= À 
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dans le cas du bloc partiellement équilibré À;,.=À, ;iet s étant uème 
associées, 


Nous garderons la notation À;, pour traiter les deux cas en mé- 
me temps. 
Estimation ‘'intra-bloc" 


Nous supposerons de façon générale, que le plan d'expérience est 
telqu'il n'y a pas d'interaction entre les blocs et les variétés. Si nous 
appelons Yi; le résultat observé en appliquant la variété i au bloc j, 
nous pouvons écrire : 


VA EN HV DU EE 


ij 
de telle sorte que y; est la somme : 


- d'un facteur ui contribution commune. 


- d'un facteur v; contribution de la variété i. 


d'un facteur b contribution du bloc j. 


- d'un facteur e; 
relation avec les autres e; 


j erreur normale, de moyenne nulle sans 


; et de variance Ce 


Pour obtenir la meilleure estimation sans distorsion de Y:; 
nous savons en appliquant le théorème de Markoff, qu'ilest nécessaire 
et suffisant de rendre minimum l'expression : 


ZT FAR NA b.) 
ij 


Différentions par rapport à vi : 


Somme étendue à l'ensemble des blocs contenant la variété i. 


expression qui s'écrit : 
- y, "2 n;;(p+ V.) +2 n,,b, 


employons la notation : 
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Différentions par rapport à b : : 


DO CORRE 0 


Somme étendue à toutes les variétés du jème bloc. 
DE D pi; =D m;(B+b;) +2 nv 


employons la notation : 


b) Y =k(L RD) ED av 


remarquons qu'enfaisant la somme des v équations a), ou des b équa- 
tions b), nous obtenons : 


CSS Eve cbr Sp, 
in j l 


L'équation b) peut s'écrire : 


k(È + 


D» 
1 
HG 
1 
M 
Le 
s« 


en rapportant dans l'équation a) : 
KY, = krv, ee HYe D nf, v. "2 (Env) 
J (l S£l 


que nous écrivons, en tenant compte des relations entre les n;;: 


s#i 


(I) (kr-r)v, -XJA,.v. =kY, - ŸJan,,y., 
j 


=kV - T=c, 


! 


où V; est le total des observations pour la variété i ; 


et T; est le total des observations pour les blocs contenant la variété i. 


CI 


Si nous appelons C le vecteur colonne : C = 
C, 


et V le vecteur colonne : V=| * 
V, 
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J 


A la matrice (a;;) avec AAA, 
a;, = r(k-1) 


Le système (I) s'écrit : 


Estimation ‘'inter-bloc. 


267 


Pour comparer les variétés en utilisant les résultats fournis par 


tout un bloc, nous écrirons comme précédemment : 


où b; est une erreur normale, de moyenne nulle, sans relation avec 


les autres b , et de variance eu ( o° + ko). 


Pour obtenir la meilleure estimation sans distorsion de Y, nous 


rendons minimum l'expression : 
UY La Le 72 2 
Dies = ky -> ni; Vi) 
j 

Différentions par rapport à 


D CR kd-ÿn;;v;)=0 
j î 


Somme étendue à l'ensemble des blocs contenant la variété i. 


ou encore : 


+ \ Ne Es 
rku > œ n;,n.) v "2 n,, Ÿ, 
comme la moyenne des Ÿ, est nulle : 


Y .=r 9% +bk =bkfi 


et, tenant compte des relations entre les n;;: 


soit : 
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où T, est letotal des observations pour les blocs contenant la variété i. 


Si nous appelons C le vecteur colonne : C = 
C, 
A la matrice À = (à,,) avec ai; = À;; 
a; mer 
Le système (II) s'écrit : 
A: Vic 
Estimation ‘'intra-inter-bloc''. 
| 1 1 
: : se PR nés dre : 
Si nous connaissons W et W O+Kko nous obtiendrons la 


meilleure estimation sans distorsion des variétés en combinant les 
estimations pondérées précédentes, Nous la noterons v;, et nous l'ob- 
tiendrons en rendant minimum : 


W 
VOUS ES DE RE LEPn UN 
ij j i 


En combinant les deux systèmes matriciels déjà obtenus nous obtenons 
le nouveau système matriciel : 


W(GRAV)+W.(C. ANI)=0 
ou-encore : 

W C + W'È = (W A+ W' À) v' 
Posons : 


P=WC+W È=Wkv-T)+W'(T-Efe u) 


où V, T sont déjà définis, et U est le vecteur colonne unité, 


Ce système s'écrit : 


(III) P =B *Y 


avec B=(b;;), b; = Wa +W' à; 
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si nous posons : 


È = r(W(k-1) + W!) 


B = W' - W 
b;;= BA;, 
b;;= a 


Lorsque nous connaissons les À;; nous pouvons grâce au système (III) 
calculer.les v.. 


Cas du bloc incomplet équilibré. 


Si nous appliquons le résultat précédent, dans le cas de l'esti- 
mation intra-bloc par exemple, au Bloc Incomplet Equilibré, les équa- 
tions (III) s'écrivent : 


ou encore en tenant compte du fait que Ÿ % = 0 


C + (nc AdE 


La différence des vŸ, est donc estimée par la différence : 


Ci - C CRC 


r(k-1)+ À TRE 


: 2k 
et la variance est : rs oc. 


Dans le cas de l'estimation inter-bloc, W=0, W'=1, 
Les équations (III) s'écrivent : 
A) 


DY, 


La différence de deux Ÿ est estimée par : 


k 
= (on PKo } 


la variance étant : 


Nous retrouvons ainsi brièvement les résultats connus de l'analyse de 
la variance du Bloc Incomplet Equilibré, à savoir que la variance est 
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la même pour n'importe quel couple de variétés. Ceci provient du fait 
que À est une constante, nous allons vérifier, dans le cas du Bloc In- 
complet Partiellement Equilibré l'influence du facteur À;.. 


Cas du Bloc Incomplet Partiellement Equilibré à m classes associées. 
rl Les À;, Un partie d'un groupe de valeurs À,, À,, ..., An avec 
> SE 


E 2 m° 


Nous allons chercher un système matriciel utilisant les paramè- 
tres de deuxième espèce. 


Faisons choix d'une variété i, et appelons (S;,) l'ensemble des 
variétés uèMES associées de i ; notons G:, la somme des vi pour tout 
JesS,,. Le système (III) s'écrit : 


PF = vf Su 


} 
u=] 


LC 


Faisons la somme de ces équations pour les n variétés qui sont uèmes 
associées de la variété i : 


zP 


iu 


sa G;,+$pA (P,: Br re. bre) 


313 


+ BA, (n, v; + Dé Gi, ae Êb, Gin) 


KUd 0 


+ BA, (Pie Gi LA TS PU RACE 


Finalement, nous obtenons les m+1 équations : 


P =a v +BVA;G; 


J=1 


ao) 
1 


CR A nEv RES É & +B ar, [a (j=1;..., m) 
k u 
où ë. est le symbole de Kroneker. 


Si nous imposons la condition : 


nous pouvons réduire ce système à un système de m équations à m in- 
connues : 
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. m-]1 : ml 
F ra ; s> BY À, G:; 1 BA, (4 SR G;,) 
u=] 


Î 
j=1 


PR Ar BD CAC 


JET 


: a : ee 
EP; =BA;n;v; 2 (aë; +8 2 AP.) G;, = 2 À, p,, (v, > G:,) 
= u= u= =] 


m-]l 
Pan cp AD ve ES [es REA DEE m| Gix 
u u=1 u 


RES RESTES oo tan le 


Système que nous écrirons sous forme matricielle : 
(IV) IMaMeeZz 


avec pour li le vecteur colonne 


P 
Il = 
2F;, » de, "Hi L 
Z le vecteur colonne 
VF 
22 — 
ïj 
et la matrice M : 
a — BA BON = AS) D roms POSER) 
M=|B(Ain,- >: Aup°) a+ f D: À (be Pi) — pu À, (pis po) 


ÉD ES AR, pr) —ctPi (pr -p".) 


m m-ilu % 
Appelons À le déterminant de M. 
Ce déterminantn'est pas nul ; le système (IV) est, en effet équi- 


valent au système (III), et | B | n'est pas nul. Si nous résolvons ce sys- 
tème de Cramer seule l'inconnue ve nous intéresse. 
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Appelons A le mineur du terme a; de la ième colonne et de la jème 
ligne, dans 


Pour obtenir Ÿ,, il suffit de faire W' = 0, W = 1 dans ve et pour V, : 
ENT), AO (ES 


L'expression de v. a été donnée par Rao C.R., dans les cas 
m=?ret-m=3. 


Estimation de la différence de deux variétés. 


à 


Nous allons traiter en même temps le cas de l'estimation ‘'intra- 
bloc'' et de l'estimation ‘'inter-bloc'', il suffira de remplacer Wpar 1 
et W!' par 0, dans le premier cas, W par 0 et W' par 1, dans le second 
cas. 


Ne Su 
À  j=0 


j+1 
ñ 


* hrs 
V: VS 


(Z He 2P,.) 
expression que nous pourrions écrire en fonction de V, et VW, si les 
variétés i et k sont tèmes associées (t{#m}), V. apparait dans : 


PROS 


i io 


ét 2 Po cars 


(nl 


P, KV, - T;, ZPx=Kk(EV-ZT) 
de même V, apparait dans P = LP, et Z P::. 


La variance de la différence est donc : 


La variance de la différence de deux variétés dépend de leur classe 
d'association, (de la valeur algébrique de A  ). 


Si ces deux variétés sont mÈèME associées, V, et V, n'apparais- 
sent que dans P. et P, respectivement, la variance de la différence est 


donc : 
2e E c: | 
A 
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Finalement la moyenne des variances des différences des deux variétés 
gécrit: 


Coefficient d'efficacité. 


Le coefficient d'efficacité d'un bloc incomplet est défini comme : 


Variance moyenne des différences de variétés dans le bloc complet 


4 ROUTE FEU DÉS TM PAT Tree 
Variance moyenne des différences de variétés dans l'estimation 
intra-bloc du bloc incomplet 


Variance intra-bloc dans le bloc incomplet = 02 


——————————————————.—.—.—— 


Variance intra-bloc dans le bloc complet = 0° 


Chaque variété étant répétée r fois, lavariance moyenne des différen- 
ces de variétés dans le bloc complet, est : 


2 
2 0: 


r 


À - Cas du Bloc Equilibré, 


ER DS a 


ee lt) 
A UMMErT Kr  k(v-1) 


E = 


B - Cas du Bloc Incomplet Partiellement Equilibré. 


Nous remplacerons W par 1et W' par 0. 


oc? 20°(v-1) A 
E = TR Do UD Lit] 
O3 2ko°|(v-1) AS nn À | 
u=] 
Ver A 
Fa 


m=l 5 
1 j+1 
HrI) ARENA 
j=1 
Remarquons que nous pouvons également définir un coefficient d'effi- 


cacité à l'intérieur de chaque classe d'association : si les variétés sont 
uèmes associées, la variance de la différence est : 
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Æ Ça - A] 0? (um) 


et le coefficient d'efficacité correspondant : 
1 À 

6) DRRON PAS AU 
rk À; A, 


Si les deux variétés sont mèmes associées, la variance de la diffé- 
rence est : 


LITE a 
AFSOS, En — 
TUE A 


Applications numériques. 


Les formules données pour v;, E sont calculables par des ma- 
chines électroniques. Lorsque le nombre de classes d'association est 
inférieur à 10, les déterminants A sont calculés en moins de vingt 
secondes. 


Pour utiliser pratiquement, un plan d'expérience du type bloc 
incomplet partiellement équilibré à m classes associées, il faudra 
avoir toutes les caractéristiques qui se sont révélées utiles dans l'ana- 
lyse, à savoir : 

V;1D,K, ren A, FE, 


u? 


Nous allons donner deux exemples. 


Exemple 1. 


Considérons le bloc incomplet partiellement équilibré à trois 
classes associées, construit à l'aide de la géométrie projective finie 
GPI(5702). 


Ve & À bu? r =6 k = 4 
ny € À n°; = 3 n, = 3 
À,=3 À, = 2 À,=3 
ON O0 OCR OS 0 
0 NUS pin 0, 200 EL ONU 
OS O0 OO 2 Où 2 


Nous allons calculer les Aj : 
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18+ 3 0 1 
A= | -1 18 - 2 -1 = 21A - À+3 4 
3 0 18+5 
se 
Af= 16 5291 M =0, A, = - 16 
23 P, CPP 
V.. = ——- 
480 
a ; 
rime E = 0,853 


Exemple 2. 


Considérons le bloc incomplet partiellement équilibré à trois 
classes associées, ayant pour paramètres : 


v = 12, b= 12, r=3, k=3 
h, =2, n #9, n, = 6, 
À,= 0, À, = 0, À3= 1. 
LOTO LAS DL 
DM O0 03 DEAR N0 Be MANOIR 
OMS 110 204 M  ) 
donc : 7 jt 
AS = 7,55 - 1.9 - 2,8 = 360 
DRE TE 
avec : A, = 55, A = 4.9, AE= 8. 
OP SP Dir 
| 360 
Nu 360 110 
as CRE NE ETS TT: 6 EC 


Analyse de certains Blocs Incomplets Partiellement Equilibrés quand 


un bloc manque, 


Comme nous l'avons expliqué dans l'introduction, le bloc est 
très souvent une unité expérimentale ; par exemple l'ensemble des ani- 
maux d'une même portée. Il peut donc arriver assez fréquement que 
les résultats de tout un bloc soient perdus. 
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Ce problème a été traité, pour les carrés latins par Yates et 
Hale, par Cornishpour un Bloc Incomplet Equilibré, par Marvin Zeler 
pour un Bloc Incomplet Partiellement Equilibré à deux classes asso- 
ciées seulement. 


Nous allons donner une solution générale. 


A - EQUATIONS DE L'ANALYSE ‘'INTRA-BLOC'" DE LA VA- 
RIANCE - 


Nous supposerons que le bloc manquant contenait les variétés vw, , 
Vo 000 M 


Nous emploierons les mêmes notations que dans le paragraphe 
précédent. Pour établir l'analyse de la variance, nous distinguerons 
deux cas: celui où la variété envisagée apparait dans le bloc manquant 
(i$ k) celui où elle n'apparait pas dans le bloc manquant, (i > k). 


Pour i>k, les équations de l'analyse "intra-bloc, sont les 
mêmes que dans le chapitre précédent : 


Dec ERIC MC 
u=]l 


u=] 


compte tenu de la relation : v; + Ÿ G, = 0 
u=]l 


Pour i<k, il faut ajouter aux résultats obtenus les résultats 
manquants. Par définition €; = kV, - T,, où V. est le total des obser- 
vations, pour la variété i, et T; le total des observations pour les blocs 
contenant la variété i. 


Il nous faudra tenir compte des observations qui n'ont pas été 
faites dans le bloc manquant, pour exprimer V, et T;. Quand la variété 
i apparait dans le bloc manquant, l'équation ‘'intra-bloc'' s'écrit : 


k m 
i + (KŸ; > Ÿ.) ES r(k- 1) Ÿ: e > À, Giu i < k 
S=] u=] 


soit en tenant compte de la relation LŸ, = 0 


a Los k Loi Lol nez 
CARVIN) SL fr (ke I) EENEINERS SANG 
s=1 j#d 
comme dans le paragraphe précédent, nous allons faire la somme de 


ces équations pour les n; variétés qui sont jèmes associées de la va- 
riété i. 
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Nous appellerons m;;, le nombre des jèmes associées de la va- 
été i, qui apparaissent dans le bloc manquant, 


Li k # LS # ES # 
Z (kv, 5 WF) EC k S; (Ke m; 2 v. 


et, en notant S, (v.) la somme des résultats pour les jèmes associées 
de la variété i, qui apparaissent dans le bloc manquant, 


2 @ V, - D A G,)=-Nnÿ $2 | «5-5 Ari, |, 
a + u 


avec x = r(k-1). Soit en tenant compte de la relation : 


7 LEE 
G;, HA SE > Gi 
u-] 
m 3 - 1 m 
-(Agn;-7 Aupi)v; +Y [aë;-ZEA (p,-p:)1G, 
u t=1 


Finalement, pour i>k, nous obtenons les m équations : 


m-1 
ee ROMA NT D VS A) 


u=] 


" 4, m-1 
= m 4 + 4 + m 
Ce el A D | «6; rs ar) | 
=1 u 


Equations que nous écrirons sous forme matricielle : 


1 1, 2E 


avec ©, le vecteur colonne %- 


8, le vecteur colonne & = 


M 1a matrice M où W' = 0 et W = 1. 


De même pour i< k: 
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3 k 
Che 
S=]l 
avec: F, = x « n 
ANNEE 11e Al, nel 
-k 0 . 0 
OMORTE 0 
et: Nu = ICE 
2 l 
0 oO 0 


En résolvant ces deux systèmes comme dans le chapitre précédent, 
nous obtenons pour V; : 


Mal N " : Fr 
AC, +9 AT (G;+kS;(vi)-m;; D v) | 


j=21 


où À, est le déterminant de ; 
et AÎ le mineur du jème terme de la première colonne. 


À, sera me déterminant de, ; remarquons que les mineurs A; sont 
les mêmes pour À, et A. 


Pour i <k, nous aurons donc : 


(Er +kS) (M) mp 5 | 


k 
S; (v;) et 2 Ve sont 
S= 


Ces équations sont encore inutilisables, car M 


j » 
inconnus. 


B - CAS OÙ LES VARIETES DU BLOC MANQUANT SONT 
mèmes ASSOCIEES LES UNES DES AUTRES - 


Nous pouvons supposer que toutes les variétés du bloc manquant 
sont mÈmes associées les unes des autres, la classe m d'association 


étant, jusqu'à présent, une classe quelconque. 


Soit une variété du bloc manquant, (i=1, ..., k) 


À 
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SM) SV -V, m,, =k-1 


1) 06 jfm, m,;,=0 


puisque les k variétés du bloc manquant sont nèmes associées les unes 
des autres. 


Ecrivons la valeur de Ÿ 


A 


PVR A (Ce Et)+ DA (Sch-miv.) 


1 2h 


AS TE = 
ACHOP ADN TG, PAZ 
j=1 
Faisons la somme de ces équations pour i=1, ..., k. 


a = LS ges j+ \ 
DAV EN PEIV = DCR in) ec; - KA va 


CE) 2 £ 
J=1 


comme À,+ kA, LE 


A 1 NAS ee j+ = 
Div. re Le Dr 0 ice | 
s A; s j=1 s 


et, en remplaçant dans l'expression de Ÿ, : 


= A DE, + = a ra \ 
ie Lfac+S are, |- La AC D DeA 0 £c;, 
à; j=1 A, A, s j=1 s 


Nous obtenons, ainsi les estimations des variétés qui apparais- 
sent dans le bloc manquant. 


Pour les variétés, n'apparaissant pas dans le bloc manquant 
(i > k) nous remplacerons les quantités S, (V;,) et Z v,, par leurs esti- 
mations numériques. 


Pour comparer la différence de deux variétés, nous devrons non 
seulement tenir compte de la classe d'association, mais encore consi- 
dérer les cas où les variétés apparaissent ensemble dans le bloc man- 
quant, où l'une des deux y apparaît, où aucune n'y apparaît. 


1/ - Nous supposerons tout d'abord que les variétés i et t sont 


mèmes associées. 


a) v, et v. apparaissent toutes les deux dans le bloc man- 
quant 
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V, apparait dans c, et ce. 


V, apparait dans €, et C;, 


donc 


varience (ŸV;, - V,) = 


b) v, et v, n'apparaissent ni l'une ni l'autre dans le bloc 
manquant : 


V: _— | 4 Gé Ÿ AY 2c;,-2 C+; +kS; (v:)- kS; RE (m;,-m,;)2 | 
1 


IE 


comme i ett sont mÈMES associées V, n'apparaissent que dans c,, et 
Vdans 


vi 
2Aï 


variance (v, = Ÿ, ) = 
ÿ à 


c) v; apparaît dans le bloc manquant, v, n'apparaît pas dans 
le bloc manquant. 


= 1 D de io A: en je 
VEN Ge A Cc;+) AÙ Zc;, -—2(A;LC, + Ÿ A! D ace) 
À, j=1 s 
1 ST 
he AC NRA (2cx; +kS;(v,)-m,,Ev,) 


1 J=1 


V, apparaît dans C. ; Ce 2 


V. apparaît dans c, ; 
. , AE 0 ? A x 
variance (v; D do te - +—EIm,AS 0° 
À, A; A, à, ds 


2/-Les variétés v, et v, sont uèmes associées, u=l,...,m-i, 
22 


a) niv, , ni v, n'apparaissent dans le bloc manquant. 


Nous supposerons que les deux variétés i et t ont M; jèmes asso- 
ciées communes dans le bloc manquant. 


me 1 x » Let , 1 £ 
V. == — EC ne) AT (Zi; - LérjtkS,(v, )-kS(v,)-(m;,-m,,)2v, 
Jj=1 
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V, apparaît dans €,, ê 


& €t dans DAA car : 


m1 : 

4 = de \ A Gi CEA A 

A, EV, À, 2 CRD À > 2 c;. 
S = S 


comme v, a m,, pèmes associées dans le bloc manquant, V, apparaît 
Smic. 


de même V, apparaît dans €,, 6;, et dans les m,, ê 


ps® 
2 o2 & —. fl LES DOS: È a Sn 
À SE Hunt : 
variance (Ÿ,-v,) = Rite > Do ietn me 
1 1 


p=1l j=1 
2 M: ; Mt) M, (m;; d | 


b) v, apparaît dans le bloc manquant, v, n'apparaît pas dans 
le bloc manquant. 


À A & m1 + A 2 mL j+ à 
D - D [as+S aire, in ne) AD S 20, | 
j=1 1 j=1 s 


ml 
1 p+l vw 
, ee RTC > A > 20) 
A — 
an A 
Mhapparalitidans 6, €, €, 5,(v.). 
Vs apparait dans €,, €, €... 
Re. Re ; cn A? 
variance (Ÿ, -Ÿ,) = (—+—) (A,- A) 02-092 - 1 52 
à, À; A, À, à; 
À À 
+— S m, Allo? + TA 'm.;, 0° 
D 119? 


o? s j+1 np+l 
Le DOC m,; m4, 
TASRRE 


ou encore : 


L 
: il 1 _. A  —. 
variance (Ÿ, - Ÿ,) = o? | (—+ —) (A°- Fi PR > A m:;) 
… à; Tu 
nee 1 F % ; 
en +—2+—ÙY ÿ AÏ a n,, M, 
LE 


2 Bi p 
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Exemple. 


Tous les Blocs Incomplets Partiellement Equilibrés construits à 
l'aide de carrés latins orthogonaux, contiennent des blocs où toutes les 
variétés appartiennent à la même classe d'association. 


Nous voulons construire un plan d'expérience, avec v variétés, 
v=pq, p< q, q étant un nombre premier, ou une puissance de nombre 
premier, 


Dans ces conditions, il a été démontré que l'on peut construire 
q(q- 1) carrés latins orthogonaux, de q lettres, c'est-à-dire q carrés 
tels que, sil'on superpose deux quelconque d'entre eux, chaque couple 
de lettres (ou de nombres) n'apparaît qu'une fois et une seule dans tout 
le carré. Comme p <q, il a été démontré que, si l'on supprime les 
p-q dernières lignes on obtient un bloc incomplet partiellement équili- 
bré, avec les paramètres : 


v=pq, b=q(q-1), k=p, r=q-1, n;=p-1,n,=q-1, À,= À,=0, À,=1. 


exemple numérique : p=3, q=d4. 
ve d2; DAT RES, MES, n, =6 


le schéma d'association est obtenu comme dans le premier chapitre, 
le plan d'expérience étant considéré comme quasifactoriel : 


11 12 13 14 
21 22 23 24 
31 32 33 34 


les 12 variétés étant représentées par un couple de nombres ij,. 


Le pian s'obtient à partir des 3 carrés latins orthogonaux en sup- 
primant la dernière rangée : 


17 LES ET LORS 4 SOS 
2 ZA ESS SAT 2 LR 2 ET 
3, STI? 43 211 2 LAS 


chaque colonne est considérée comme un bloc, la variété est obtenue 
en ajoutant le numéro de la rangée au chiffre figurant dans le carré, 
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COR 2e) NTI 215 1e l'A os 22) 
7 Ne A AS, 21 22)  (ld 22 151) 
2 12 20 51) 0 (1532-22 034) (140 21252) 


LD 0 O0 2 PQ A 
cbr, 21000005 GERS MOPRE EX per] ÉUG-eE 
OS 201 ln 
6+1 1 ï 
= _ = Ce Lie 
he 1 6+1 LA GO A EST À SGA 26 
2 2 6+2 
A, = 198 
6 
variance (v, -v,)> RE loee Ft 5; 


en supposant que le bloc (11 22 33) manque. 


les variétés 12 et 32 sont premières associées, elles n'apparaissent 
ni l'une ni l'autre dans le bloc manquant : 


posons i=12, t=32 


donc, en appliquant la formule trouvée précédemment : 


x 26 (ei 
variance (v, - V2) deco (60.5) PET 


EIPSE 
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NOMBRES ALÉATOIRES, SUITES ARITHMÉTIQUES, 
MÉTHODE DE MONTE-CARLO 


J. BASS 


SOMMAIRE 


La méthode de Monte-Carlo est une méthode de calcul numérique fai- 
sant appel à des analogies probabilistes. Elle utilise des nombres aléa- 
toires qui sont Les résultats d'épreuves successives sur une variable 
aléatoire donnée. Les nombres aléatoires Les plus simples sont ceux dont 
la distribution est uniforme entreoet 1. IL en existe des tables, cons- 
truites par des procédés difficiles à justifier. On montre, en prenant 
comme modèle le calcul d'une intégrale définie, qu'il ya intérêt à subs- 
tituer aux méthodes probabilistes des méthodes arithmétiques. Ces métho- 
des reposent sur les propriétés des suites de nombres réels uniformément 
denses dans un intervalle. Les théorèmes de Æ.heyl en précisent La struc- 
ture et en donnent des exemples utilisables. Le plus simple est celui 
des multiples modulo 1 d'un nombre irrationnel 4. La discussion de l'er- 


5 N-1 

reur commise en remplacant l'intésrale simple HXJAx par + sS f(Xn), 
0 nx0 

où X, = An, mod 1, montre que la convergence est plus rapide qu'avec des 

nombres aléatoires véritables. On établit une formule permettant de cal- 

culer dans ce cas une limite supérieure de l'erreur commise. Lorsque A 

est aléébrique, et moyennant des conditions de résularité sur f(x), des 


formules de Richtmyer montrent que l'erreur est proportionnelle à+.une 


discussion est amorcée pour discuter l'ordre de $randeur de l'erreur lorsque - 
X,= An? mod 1. La précision semble moins bonne, mais la suite An? pré- 
sente cependant quelques avantages sur La suite An: elle n'a pas de cor- 
rélation interne, et elle peut ètre utilisée pour former des nombres 
aléatoirés prenant les valeurs O et 1 avec des probabilités éé$ales. 


1 - INTRODUCTION - 


Donner une définition générale de ce qu'on appelle méthode de 
Monte-Carlo est une entreprise difficile. Il semble qu'il s'agisse 
d'une méthode de calcul analogique, dans laquelle on remplace un 
phénomène naturel, présentant certains caractères aléatoires ou sta- 
tistiques, par une image de nature probabiliste et mathématique de 
telle sorte que, par des opérations statistiques pures, on puisse 
simuler le phénomène physique aléatoire, et en calculer les para- 
mètres d'une façon plus économique que par des expériences directes. 
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Une telle définition est très vague, et je n'ai pas l'intention de 
la discuter ici. Je me limiterai à un point de vue plus restreint. Je 
considérerai comme ‘méthode de Monte-Carlo" toute méthode qui 
permet de résoudre numériquement un problème mathématique par 
une analogie probabiliste. Le ‘problème mathématique” consiste à 
calculer un nombre ou un ensemble de nombres définis par des opé- 
rations d'analyse pure, où n'apparaît pas la notion de probabilité : 
calculs d'intégrales simples ou multiples, inversions de matrices, 
résolution d'équations intégrales, et cela sans qu'intervienne l'origine 
physique du problème. Il s'agit de découvrir, entre le formalisme 
du problème posé, et celui d'un certain problème de calcul des pro- 
babilités, une analogie telle que le second puisse servir à résoudre 
le premier. Or, résoudre un problème de probabilités, cela se tra- 
duit en pratique par une expérimentation statistique qui fournit des 
valeurs approchées du résultat, avec une exactitude d'autant meilleure 
que le nombre d'épreuves est plus grand. D'où des méthodes pratiques 
de calcul numérique, dans la mesure où l'on dispose de moyens ma- 
tériels permettant d'effectuer, avec une précision suffisante, un 
grand nombre d'expériences successives. 


La statistique fournit des fréquences et des valeurs moyennes. 
La méthode de Monte-Carlo consiste donc à interpréter, après des 
transformations convenables, le problème posé comme un calcul de 
fréquences ou de moyennes. Mais on peut toujours traduire une fré- 
quence en langage de valeurs moyennes. D'une façon un peu schéma- 
tique, on peut donc dire que le but de la méthode de Monte-Carlo est 
le calcul des valeurs moyennes, c'est-à-dire d'intégrales simples ou 
multiples. 


Bien que cette conclusion soit loin d'être absolue, je vaisla 
choisir comme base de départ, et je prendrai dans ce qui suit com- 
me modèle de problème le calcul numérique d'une intégrale simple 
ou multiple. Si le domaine d'intégration est borné, on peut toujours 
leP supposer intérieur à l'hypércube CO NON ET ONE 
moyennant une transformation linéaire préalable. Si la fonction à 
intégrer est définie dans un domaine intérieur à cet hypercube, on 
peut toujours la prolonger à tout l'hypercube en lui donnant la valeur 
0 dans la partie de C où elle n'est pas initialement définie. On est 
donc conduit à l'intégrale 


I “ ÉÉCSVEON, POLATARE 
Are 
I - CALCUL D'UNE INTEGRALE - 


Soit f(x) une fonction bornée, intégrable, définie sur le segment 
(0. 1). L'intégrale 
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ei 
I = | f(x) dx, (210) 
«0 
qui a priori n'a aucune relation avec le calcul des probabilités, peut 
être considérée comme la valeur moyenne (espérance mathématique) 
d'une variable aléatoire X uniformément distribuée entre 0 et 1. 


D'après la loi des grands nombres, si x,,..., x, sont les ré- 
sultats de N épreuves successives sur X, la moyenne arithmétique 


1 
CT N f(x.) +... + f(x)] (2.2) 


tenden probabilité vers I lorsque N —>« , Il est d'ailleurs commode 
de considérer x,,..., x, comme les résultats d'épreuves successives 
sur N variables aléatoires indépendantes uniformément distribuées 
SUR (OS 1): 


Désignons par © l'écart-type de f(X), défini par 


L 


GÉ,= vi SR - | HOT dx - I’. (2.3) 
. «/0 


e) 


On sait que S, est une variable aléatoire qui a pour moyenne (stochas- 


(el 
\N 


ment distribuée suivant une loi normale réduite. 


tique) I, pour écart-type et que ÿ (S,- 1) est asymptotique- 


Par suite, si N est suffisamment grand, 


Probe - Dhs 2 / e 2at 
h 


V2 4/nVx 


Ho = 1! est la valeur absolue de l'erreur commise sur l'évaluation de 
I nr on prend $, comme valeur approchée. Si l'on fait 4N épreuves 
au lieu de N, on a la même probabilité de commettre sur I une erreur 
deux fois plus petite. C'est ce qu'on exprime en disant que l'erreur 


7 1 : : ; : 
décroît comme NX - Si l'on a fait une évaluation grossière de ©, et 


qu'on se fixe un seuil de probabilité, on peut calculer le nombre de 


termes N nécessaire pour obtenir I avec une précision donnée à l'a- 
vance. Si par exemple on néglige des éventualités de probabilité 0,05, 
le nombre de termes N permettant d'obtenir I à 0,001 prèsest défini 
par 


ie = 0,05 avec h = 0,001, 
= 
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soit h\N = 1,96, et approximativement N = 4.10°0°. 
o 


Mais ce résultat n'a pas grande valeur en soi, puisqu'il dépend du 
seuil de probabilité accepté, qui est tout à fait imprécis. 


1 
Le seul résultat à retenir est que l'erreur décroît comme NN > 


c'est-à-dire lentement. La suite des x, est irrégulière et remplit 
progressivement l'intervalle (0.1). Mais on ne sait rien d'autre sur 
elle que son caractère aléatoire. Elle ne possède aucune organisation 
qui, sans modifier les principes, puisse favoriser la convergence. 


Il est évidemment maladroit de calculer une intégrale par la 
formule (2.1). Les méthodes numériques usuelles sont plus efficaces. 
La ‘méthode des rectangles", moins précise que la méthode des tra- 
pèzes ou les méthodes d'interpolation par polynômes, donne déjà une 
précision en & . La méthode statistique consiste à prendre des points 
au hasard, de sorte que leur répartition, médiocre pour les petites 
valeurs de N, en ce sens qu'il peut y avoir peu de points dans un in- 
tervalle, et beaucoup dans un autre, ne s'uniformise vraiment que 
lorsque N devient grand. Au contraire, la méthode des rectangles con- 
siste, une fois N fixé, à placer N - 1 points suivant la meilleure ré- 
partition possible entre 0 et 1, c'est-à-dire de telle sorte que leurs 
abscisses soient en progression arithmétique. Une valeur de N rela- 
tivement petite peut donc déjà donner une précision que la méthode 
probabiliste n'atteindrait qu'avecune valeur bien plus grande. 


On notera cependant que l'évaluation de la précision dans la mé- 
thode des rectangles n'est possible que sila fonction f(x) est suffisam- 
ment régulière. Si f(x) admet une dérivée continue, telle que, sur 
l'intervalle (0.1), on ait |f'| < M, on démontre à l'aide d'une intégra- 
tion par parties que 


LAS 
MINES — 
TK 16) 


k=0 


€ 5N (2.5) 


Cette formule n'est donc pas applicable sif(x) n'a pas une dérivée con- 
tinue, ou est définie empiriquement. On est alors conduit à choisir N 
d'une façon plus subjective que rationnelle, et on doit évaluer la pré- 
cision du calcul par des procédés empiriques. Si la précision semble 
insuffisante, on doit recommencer l'ensemble des opérations avec une 


valeur supérieure de N. Il est d'ailleurs difficile d'évaluer à vue la 
précision. 
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Figure 1 


Au contraire, la méthode statistique ne nécessite pas un choix 
initial sur N. Au fur et à mesure de la croissance de N, la valeur de 
S, se stabilise et on la voit tendre vers la limite I. On s'arrête quand 
on a atteint la précision souhaitée. 


HI - NOMBRES ALEATOIRES - 


Le premier problème qui se pose alors est de construire une 
suite de ‘nombres aléatoires”. Il s'agit de nombres aléatoires au sens 
restreint, c'est-à-dire de nombres dont la suite ordonnée puisse être 
considérée comme le résultat d'une suite d'épreuves sur une variable 
aléatoire uniformément répartie entre 0 et 1. Je ne m'occuperai pas 
ici du problème général de la construction des nombres aléatoires ayant 
une répartition quelconque. Si X a pour fonction de répartition F(x), 
la variable aléatoire Y = F(x) est uniformément distribuée entre 0 et 
1. En effet, puisque F(x) est non décroissante, 


Prob [X < x] = Prob [Y < F(x)] = F(x), 
et si l'on pose F(x) = y, on a bien 
Prob: [Ye yi.=.y 
On peut donc toujours ramener une variable aléatoire quelconque 


à une variable uniformément distribuée entre 0 et 1 et le problème qui 
se pose est un problème purement pratique. 
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Pour construire une succession de valeurs numériques qui soient 
le résultat d'épreuves indépendantes sur une variable aléatoire unifor- 
mément répartie entre Oet 1, on doit recourir à un jeu statistique con- 
venablement choisi. On obtient ainsi ce qu'on appelle une table de 
nombres aléatoires. Il n'est d'ailleurs pas facile de réaliser un jeu 
statistique qui soit à la fois assez précis, assez rapide et assez garanti 
pour donner avec un nombre de décimales adapté aux machines à cal- 
culer modernes plusieurs dizaines de milliers de nombres successifs. 
C'est ce qui a été fait pour calculer les tables RAND, et les tests 
statistiques semblent prouver que le résultat est correct. 


Cependant d'autres méthodes ont été proposées pour calculer 
des nombres aléatoires. Ces méthodes reposent sur des calculs de 
nature arithmétique dans lesquels le hasard ne joue aucun rôle. Ces 
calculs sont organisés de telle sorte que la suite x,, x,,... qu'ils 
fournissent soit compliquée et donne l'impression du hasard. Les 
nombres obtenus sont ensuite testés par des procédés statistiques. 


Malheureusement, les méthodes arithmétiques sont souvent dif- 
ficiles à étudier du point de vue mathématique. Comme les tests sta- 
tistiques ne sont pas des preuves absolues de la validité des hypothèses, 
les méthodes arithmétiques généralement employées restent très in- 
certaines. On ne peut pas affirmer qu'elles puissent jouer le même 
rôle que les méthodes de nature spécifiquement probabilistes. 


Il existe cependant certains cas où l'arithmétique est une science 
assez avancée pour qu'on puisse l'invoquer en toute connaissance de 
cause lorsqu'il s'agit de construire des nombres aléatoires. C'est à 
l'examen de ces théories arithmétiques qu'est consacrée la suite de 
cette étude. On notera que, dès le moment où l'on utilise des théorèmes 
précis, les méthodes de calcul qu'on obtient n'ont plus rien de statis- 
tique, sauf l'apparence. Ce sont des méthodes mathématiques, pré- 
sentant de simples analogies de forme et d'emploi avec les méthodes 
statistiques. On notera aussi qu'elles présentent un avantage, c'est 
que les “rombres aléatoires" arithmétiques n'ont pas besoin d'être 
rassemblés et prélevés dans des tables. Ils se calculent par des pro- 
cédés compatibles avec la technique habituelle de programmation des 
machines à calculer. 


IV - THEOREMES DE H.WEYL(U) - 
Donnons d'abord une définition. Soit 
>. 


no... 


(1) Référence (5). 
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une suite de nombres compris entre 0 et 1. Considérons un intervalle 
(a, B)tel que 0O0<a <B <1. Supposons que, parmi les N premiers 
termes de la suite, ily en ait N' dans l'intervalle (a, B ). Si le rapport 


N' 
NN tend, lorsque N— « , vers une limite égale à B -a , et cela quel 


que soit l'intervalle (a, B ) choisi, on dit que la suite x, est uniformé- 


ment dense, ou équirépartie, sur (0.1). 


Les théorèmes de H. Weyl s'énoncent de la façon suivante : 


THEOREME 1- 


Pour que la suite x, soit équirépartie sur (0.1), il faut et il suffit 
que, pour toute fonction f(x) intégrable au sens de Riemann sur (0.1), 
on ait 
er 


| en 
0 N > œ@ N n=<0 


THEOREME 2 - 


Pour que la suite x, soit équirépartie sur (0,1), il faut et il suffit 


que, pour tout entier 1 f 0, on ait 


s e2 7% = 0 


Rs 

pu Nr 

On définit d'une façon analogue une suite de points 
EU SR SN RS 


équirépartie dans l'hypercube fondamental C de l'espace à dimensions, 
c'est-à-dire l'hypercube dont les points P ont des coordonnées toutes 
comprises entre 0 et 1. Soit ® un domaine quarrable intérieur à C, et 
de volume V. Supposons que, parmi les N premiers points P,, ily en 
! 

ait N' à l'intérieur de ®. Si le rapport — tend lorsque N—> « vers 
une limite égale à V, et cela quel que soit le domaine ®, on dit que la 
suite P est équirépartie dans C. On démontre alors que : 


THEOREME 3 - 


Pour que la suite P_ soit équirépartie dans l'hypercube C, il faut 
et il suffit que, pour toute fonction f(P) intégrable au sens de Riemann 


dans l'hypercube _C, on ait 


296 J. BASS 


où dw est l'élément de volume de l'espace. 


Les théorèmes 1et3 permettent donc de calculer numériquement 
des intégrales simples et multiples. Le théorème 1 est l'équivalent 
arithmétique exact du théorème statistique du paragraphe 2. La seule 
différence, c'est que, en arithmétique, il y a convergence au sens or- 
dinaire du second membre vers le premier, alors qu'en statistique il 
y à convergence en probabilité (ou convergence presque sûre), ce qui 
est moins précis. 


La raison "physique" de validité de ces théorèmes est la même 
dans les deux cas. La suite x, tend à remplir ‘d'une façon homogène" 
le segment (0,1). La suite P, tend à remplir "régulièrement" l'hyper- 
cube C. La moyenne des valeurs de f en des points régulièrement ré- 
partis est une approximation de la valeur moyenne théorique, c'est-à- 
dire de l'intégrale. Le sens du mot “remplir” ne doit naturellement 
pas être interprété d'une façon trop concrète. La suite x, est dense 
dans le segment (0.1). Il existe des points de la suite x, aussi voisins 
qu'on le désire de tout point entre 0 et 1. Mais cette suite est dénom- 
brable et de mesure nulle. Elle ne recouvre pas (au sens mathémati- 
que du mot) le segment (0.1). Cependant, si l'on considère une suite 
de segments de longueur € ayant pour centres les points x,, elle re- 
couvre effectivement le segment (0.1), et cela si petit que soite . 


Le problème du calcul numérique des intégrales par la métho- 
de de Monte-Carlo est donc transformé en un problème arithmétique. 
Il s'agit maintenant de former des suites équiréparties, et d'étudier 
l'approximation qu'on doit attendre en remplaçant 


? 


dl f(Pjduw par & S f(E) 


NzQ 


V - SUITES EQUIREPARTIES - ETUDE DE LA SUITE An - 


Le théorème 2 du paragraphe 4 permet de démontrer facilement 
que . 


THEOREME 4 - 


Si À est un nombre irrationnel, la suite An des multiples de 2 
est équirépartie modulo 1, 


Cela veut dire que la suite x, des parties fractionnaires de An est 
équirépartie sur (0.1). 


En effet, 
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(522) 


nzQ 


Si 1 est un entier non nul, 1A est irrationnel et e? 7'"- 1 n'est pas nul. 
Le numérateur est majoré par 2. Le second membre tend donc vers 0 

1 
avec N- 

Pour représenter la suite 
An modulo 1, il est commode 
de considérer 2rAn comme 
l'abscisse curviligne d'un point 
sur une circonférence de rayon 
1. Si À est irrationnel, on dé- 
montre facilement que les 
points ayant pour abscisses 
curvilignes 2rA, 4nA, 6rA,... 
sont denses sur le cercle. Ce 
qui est remarquable et moins 
élémentaire, c'est qu'ils le 
soient uniformément, au sens 
indiqué ci-dessus. 


Figure 2 La somme 


ce e? ir lan 
À 
se construit d'une façon analogue. 
C'estune somme de nombres com- 
plexes de module 1 qui, dans le plan 
complexe, sont représentés par des 
vecteurs unitaires tels que chacun 
fasse avec le précédent un angle 
constant égalà 2n1A. Ces vecteurs 
sont les côtés d'un polygône régu- 
lier non fermé, inscrit dans un cer- 
cle de rayon ne se 
Y 2 sinrlA 
On notera que An ne s'annule 
jamais mod. 1, sauf pour n=0, 
mais qu'il existe des valeurs den 
pour lesquelles An est aussi proche 
qu'on le désire de 0, mod.]1,. 


Figure 3 
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Si, au lieu d'un nombre irrationnel À, on introduit une fraction 
irréductible = , le théorème de H.Weyl n'est naturellement plus exact. 


Nous verrons plus loin dans quelle mesureil est approximativement 
valable, 


La suite 


est périodique de période q. Car 


2ime (n+q) 2imên 2iTp 2im£n 
q q = q 


e MRC e ne 


Elle n'est plus dense sur le cercle, car elle n'a qu'un nombre 
fini de points, sommets d'un polygône régulier fermé. Ces points sont 
tous distincts. Ceux qui sont le plus près de l'origine sur le cercle 
trigonométrique correspondent à des valeurs de n telles que pn = +1, 


2ir 
mod.q. Ils ont pour images les points e 4 , Pour trouver la valeur 
de n telle que, par exemple 


pn=1i mod q 


on utilise l'algorithme 
d'Euclide. On fait la 
suite de divisions qui 
permet de trouver le plus 
18/8 n=6) grand comme diviseur de 


(m=1) ii 
pet q, lequel est ici égal 
3/8 9/8 (n=3) ail. . 
Exemple. 
12/8 
(n=4) Etude de la suite 
8 
31. On a 
L 15/8 (n=5) 
à 8=2x34+2 
6/8 (n = 2) 3=22%x1+1 
d'où 
Figure 4 8-3=2x3-1 
8=3x3-1 


Sin = 3, on a bien 


3n = 1 + mult 8 
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VI - SUITES DE H. WEYL - CORRELATION DANS LES SUITES 
ARITHMETIQUES - 


De la suite An, on peut déduire par récurrence d'autres suites 
équiréparties. H. Weyl a démontré le théorème suivant : 


THEOREME 5 - 


Si p(t) = At’+A;t”7 +... + A, est un polynôme de degré v>1, 


dont au moins un coefficient, autre que À,, est irrationnel, la suite 
p(n) est équirépartie modulo 1. 


On en trouvera dans [1] une démonstration utilisant la notion de 
fonction de corrélation. 


Ce théorème met donc à la disposition du calculateur un grand 
nombre de suites qui peuvent jouer le rôle de suites de nombres aléa- 
toires. Le problème qui se pose est de savoir quelles sont parmi ces 
suites les plus favorables au calcul. Quelle est l'influence du degré v 
du polynôme et de la nature des nombres irrationneis ? Quelle est 
d'autre part l'influence dufait que, dans la pratique, on remplace obli- 
gatoirement tout nombre irrationnel par une approximation rationnelle ? 


Avant de répondre à ces questions, nous allons examiner un peu 
plus en détail la nature des suites équiréparties qui sont fournies par 


les théorèmes de H. Weyl. 


Nous aurons besoin de nous servir d'un lemme, qui est dû à H. 
Weyl, et qui constitue le 


THEOREME 6 - 


Pour que la suite de points P,, ayant pour COOrdOnnÉées x, y ,12,,..- 
soit équirépartie dans l'hypercube C, il faut et il suffit que la suite 


Mere over Lz, EME. 


soit équirépartie mod. 1 sur le segment (0. 1), quels que soient les en- 
tiers 1,, 1,, 1,,..., non tous nuls. 


Nous allons déduire du théorème 6 que 
THEOREME 7 - 
Si x, = p(n) mod 1 est une suite satisfaisant aux hypothèses du 


théorème 5, avec v > 2, la corrélation entre x, et x,,, est nulle pour 
tout hentier > 1. 
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Démonstration. 


Montrons d'abord que la suite x, Xnn est équirépartie dans le 
carré C. D'après le théorème 6, il suffit de vérifier que la suite 


1 1Xn 1" 1x,» 


est équirépartie mod.1 pour tout couple d'entiers 1,, 1, non simulta- 
nément nuls. Cela résulte de ce que 1,@{t) + 1,p(t + h) est un polynôme 
de degré v en général, de degré V- 1 si li +12 = 0, satisfaisant aux 
hypothèses du théorème 5. 


D'après le théorème 3, on a alors 
N-1 


} 1 à Ë ®,» + 
NU [/ xy dx dy à 


ont pour moyenne _ la moyenne du produit 


6-2)6-3) 


estnulle. Le coefficient de corrélation (envisagé au sens des moyennes 
arithmétiques) entre x, et x,,, est nul. Nous dirons que la suite x, ne 
présente pas de corrélation interne. 


Comme x,et x 


n+h 


Remarque. 
Ce résultatn'est pas exact dans le cas de la suite An. Désignons 
par x, la partie fractionnaire de An, par y, celle de A(n +h), par R 
celle de Ah, On a 
VA EC TE SEC TER 
VE GRR Et L'S TTRRR PERR RTE 
Y, est donc relié à x, par une relation indépendante de n. Le point 


(x, y )se trouve sur l'un des deux segments dessinés sur la figure 5. 
La suite (x,, y,) n'est pas dense dans le carré C. 


On a d'ailleurs 


: NS Ù Ÿ ( R) L >> 
N À * - FD ER XX =) 
\ LE De n <IR e. 
os 1 
= IN SENS 
N ul x,(x, + R) N <= X,,. 


n=0 
Xh? 1=R 


NOMBRES ALÉATOIRES, SUITES ARITHMÉTIQUES 301 


La suite x, étant équirépariie, 
la première somme est égale à 


F. PUR 
+ = =+— 
di xX(x + R)dx RhAICSÉ 


La seconde est égale à l'inté- 


1 
DA grale de la fonction f(x) définie 
R par 


f(x) = 0 si x <1-R 


IX STOXX CTI 


Elle vaut 
X 
0 OR 1} »1 FR 
1 CPE Cu 
«/ 1=R é 
Figure 5 X\, ÿYh à donc pour moyenne 
2 
TERRES 
SZ 2 


1 ; 1 
Comme x,et y, ont pour moyenne 2° pour moyenne quadratique 3” le 


coefficient de corrélation entre x, et y, a pour valeur 


3 
DR) 


1 
Il est compris entre - 8 et _ et ne s'annule que pour des valeurs 


exceptionnelles de h. 
VII - SUITES A DEUX VALEURS - 


En vue des applications numériques, (problèmes de cheminements 
discontinus), il est important de connaître des suites arithmétiques à 
deux valeurs, jouant le rôle de suites de nombres aléatoires à deux 
valeurs également probables. Ces suites ont autant d'applications que 
les suites équiréparties. 


En arithmétique comme en statistique, les suites équiréparties 
permettent de former des suites à deux valeurs. Enstatistique, si X 
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estune variable aléatoire équirépartie sur (0.1), la valeur aléatoire Y 


égale à 0si0<X< " Et Si 5 <X <1, prend bien deux valeurs 


2 : 
équiprobables 


En arithmétique, on démontre le théorème suivant : 
THEOREME 8 - 


Soit p(t) un polynôme de degré v > 2 ayant un coefficient irra- 


tionnel (autre que le terme constant). La suite des parties entières de 
(n) prend modulo 2 les valeurs 0 et 1 avec les fréquences 1/2 et 
=] Ex Elle n'a pas de corrélation interne. 


: À 1 
Désignons en effet par x,, y, les parties fractionnaires de 2 p(n) 


et de 5 p(n+h), où h est un entier positif. 


On peut démontrer (cf. [1], p.26) que la suite des points (x,, y.) 
est équirépartie dans le carré C du plan. Par homothétie, on en déduit 
que la suite (2x,, 2y,) est équirépartie dans le carré C' défini par 


O7 0 ve< 2. 


Désignons par &,, n, les parties entières de 2x,, 27%, prises 
modulo 2. £, est égal à 0 ou à 1 suivant que la partie entière de 2x, 
est paire ou impaire. Le point (£,, n,) peut occuper quatre positions, 
qui sont les sommets d'un carré et ont pour coordonnées 


(0. 0), (01); (1.0), (LES) 


Les droites x = 1 et y = 1 séparent 
le carré C' (figure 6) en quatre car- 
rés C , Cis Cros Ce 


Lorsque le point 2x,, 2y, est par 
exemple dans le carré C,, le point 
(E., n,) occupe la position (1.0). Or 
la probabilité"! (au sens de fréquence 
arithmétique) pour que le point (2x,, 
2yn)soit dans C, est 1/4. C'est aussi 
la "probabilité" pour que le point 
Figure 6 (Ë,, n,)se trouveen (1.0). Le même 
raisonnement s'appliquant aux trois 
autres positions, on voit que les 
quatre positions du point (E,, n,) ont 
la même fréquence 1/4. On voit 
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d'autre part sans difficulté que E,, partie entière de 2 x,, est aussi la 
partie entière de p(n). 


Le produit E,n,=E,&,,, prend, avec des fréquences égales, les 
quatre valeurs 0, 0, 0, 1. Sa moyenne vaut donc 1/4. Comme E, etn,. 
ont pour moyenne 1/2, le coefficient de corrélation entre £, et E,,, est 
nul, ce qui démontre le théorème. 


VIII - APPROXIMATIONS DANS LE CALCUL D'UNE INTEGRALE 
SIMPLE - 


Nous nous proposons ici de chercher une limite supérieure de 
l'erreur commise en remplaçant l'intégrale 


/ Med (8.1) 
e 0 


par 


N-] 
S fx,) (8.2) 
n=0 


AE 


où x, est la partie fractivnnaire de An, et où À est un nombre positif 
donné. Nous savons que l'erreur tend vers 0 lorsque À est irrationnel. 
Cependant, dans la pratique, on doit naturellement remplacer le nom- 


: $ : ; . p > 
bre irrationnel par une fraction irréductible : . Nous examinerons 


P 
d'abord le cas théorique À irrationnel et ensuite le cas pratique À = : : 


Nous aurons besoin dans la démonstration de résultats élémen- 
taires concernant les fractions continues(1), Soit A un nombre positif 
non entier, À sa partie entière. On pose 


n 1 
AV= A +— a, > | 
a 
di 
a = à D > il 
1 1 a, 4 
sat 1 
a, = à, 7 a, > 
3 
etc. 


On forme par ce procédé une suite a, à,, à ,... de nombres 


(1) Consulter à ce sujet: G. Valiron, Théorie des fonctions (Masson), p.17. 
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supérieurs à 1 dont les parties entières sont appelées quotients in- 
complets. Si À est rationnel, il arrive au bout d'un nombre fini d'opé- 
rations que a, devienne entier. Si À est irrationnel, a, n'est jamais 
entier, et le calcul se poursuit indéfiniment. Si l'on remplace a, par 
sa partie entière, on obtient une expression approchée de A, par des 
opérations rationnelles sur des entiers. 


m 


Qn ? 


Cette représentation approchée est une fraction irréductible 


appelée réduite. On démontre les résultats suivants : 


1/ Lorsque m——> «, P, et Q, augmentent indéfiniment. 


2/ Fn_ tend vers une limite égale à A. 


Qn 


P9 
@, 


3/ | A- 


< 


4/ La suite Q,est crois- 
sante et Q,,, > 2Q. Si Aest 


: , E 
rationnel, la suite Q. a un 
m 


KT£ (1) ei 
NN nombre fini de termes. À par- 


tir d'un certain rang, = reste 


Qn 
IN de Ceci posé, nous allons 


d'abord faire l'hypothèse que 

f(x) est une fonction continue , 
Forme primitive de f(x) et qu'elle est représentable 
entre 0 et 1 par une série de 
Fourier absolument convergen- 
te. Remarquons que ce quiest 
continu, c'est non pas seule- 
ment f(x), fonction définie sur 
(0.1), mais bien la fonction pé- 
riodique engendrée par f(x). Il 
en résulte que nécessairement 
f(0)=f(1). Cette condition n'est 
pas toujours vérifiée dans la 
pratique. Il faut alors faire 
subir à f(x) une modification 
qui la rende continue sans que 
le calcul de l'intégrale s'en 
Figure 7 trouve perturbé. On représente 


Forme modifiée de f(x) 
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par exemple f(x) à une échelle des abscisses plus petite, c'est-à-dire 
sur un intervalle (0, a), a < 1, et l'on prolonge f(x) par une fonction 
linéaire, comme l'indique la figure 7. 


On considère donc la fonction f (x) quiest égale à (&) si O<x<a, 
et qui varie linéairement de f(1) à f(0) lorsque x varie de a à 1. On a 


D ea 2, 1 es 
/ f, C9dx = / (Æ)ax+ — [f(1)-f(0)] = af f(x)ax+ = [f(1)-f(0)] 
"0 e’0 0 
Le calcul numérique de l'intégrale de f (x) fournit la valeur de l'inté- 


grale cherchée. 


On peut aussi doubler f(x) par symétrie, c'est-à-dire introduire 
(figure 8) sur l'intervalle (0.1), la fonction g(x) égale à 


f(2x) si OxEE 


f(2 - 2x) si 5 < a il 


On a 


Dir 


1] 


j 


On remarque que, sif(1) #f(0), 
ce qu'on a supposé,f(x), con- 
sidérée comme fonction pério- 
dique, n'est pas une fonction 
paire, alors que g(x)est néces- 
sairement paire. Son dévelop- 
pement en série de Fourier n'a 
que des termes en cosinus, dont 
seuls subsistent ceux dont le 
le rang est pair. 


1 » 
g(x)dx = 2] f(2x)dx = 
+/0 


- [° rooéx 


Supposons d'abord que 
f(x) soit un polynôme trigono- 
métrique ayant un nombrefini 
de termes : 


Figure 8 f(x) =D ere 


-r 


C =c, (Conjugué de c). 
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Comme J' tenax = ©, l'erreur commise dans le calcul de l'intégrale 
vaut 0 


1 es 2 imkAn 
IRON »ù ce 
N k 
Ne k4#0 
1 bai 
= es 2 imkAn 
FN CR ÿ e 
k£#0 n=0 


La somme qu'on doit diviser par N comporte un nombre fini de 


1 ; : 
termes bornés. L'erreur tend donc vers 0 comme += . Si l'on désigne 


N 


par k, celui des entiers 1, 2,..., r pour lequel e? MAS tte plus petit 
module, on a 
1 ä : 
INT <R Lez 1] 2 cul (8.3) 


k|=1 


THEOREME 9 - 


L'erreur commise en appliquant au calcul de l'intégrale d'un 
polynôme trigonométrique la méthode de Monte-Carlo, sous la forme 


arithmétique que justifie le théorème de H. Weyl, avec la suite x, = An 
mod 1 (A irrationnel) tend vers 0 comme KR ‘ 


Le résultat est donc meilleur qu'avec des nombres aléatoires 
véritables, car alors l'approximation est en , etnonen à . Mais 
en pratique le coefficient de à dans la formule (8.1) risque d'être 
grand, à cause du dénominateur. Il faut donc dans les deux cas beau- 
coup de termes (grande valeur de N) pour obtenir une bonne précision. 

Si maintenant la série qui représente f(x) est infinie, on a 


f(x) = Y cet TS (8.4) 


Le) 


On va supposer qu'il existe deux nombres positifs M et À tels que 


(8.5) 


La série est alors absolument convergente. On sait que, si f(x) est 
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continue et si f'(x) existe et est continue, sauf en un nombre fini de 
points où f'(x) a des sauts, alors À = 1. Si f' est partout continue et 
si f'' existe eta les propriétés précédemment requises pour f!', À =2, 
etc. le reste de la série |c,;| est majoré, pour k f 0, par 


es 2M 
RM 8. 
TN AT 
L La différence entrel'intégrale 
er et sa valeur approchée est égale à 
N 1 ik A 
_ et 2 irkan 
JE & CxN ee e 
1 e 2 iTTkKAN _ il 
D ul 
0 k-1 k k+i x KO 


: iTkA 
Les dénominateurs e? 7" - 1 


Figure 9 ne sont jamais nuls. Mais nous sa- 
vons que, lorsque k— «w, leur 
module prend des valeurs aussi 

voisines de 0 qu'on le désire. Pour majorer I,, séparons les termes 
de rang supérieur à Q, et les termes de rang inférieur. 


Les premiers ont une somme inférieure en module à 


2 > le, (8.8) 
kK=Q+1 
fl Ke e2 irmkaAn 
car N à 


est en module inférieur à 1. 
Sur ces termes, on commet donc une erreur inférieure à 


R, ne (8.9) 


ire . PAC iTkA 
Désignons maintenant par n le minimum de | e? ”*" - 1| lorsque 


k prendles valeurs 1, 2,...,Q. n est une fonction de Q. Pour évaluer 


à F 
n, utilisons la représentation approchée de A par des HEURE : 
On peut toujours mettre À sous la forme 


18 Ÿ 


TARN En Ton. 8.10 
MO TO RreeS Pas Ge) 
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Si l'on néglige le terme en %,, la suite des nombres complexes 
2 rh 
e 
devient périodique. Si k parcourt l'ensemble des valeurs 


ik Ppocos Q,_» 


Pm 


Om 


2 ik 
la plus petite valeur de | e = 1| est égale à 


217 


- 1 (824100) 


On sera sûr de ne pas dépasser cette valeur si l'on choisit pour 
Q un nombre inférieur ou égal à Q,, . Choisissons Q=Q, ..Comparons 
alors la suite des points 


En 
2 imk Q 
e 


m 2 irkA 


et e 


sur le cercle trigonométrique. On a 


KES k®, 
kA = Hess. 
Qn Qn Qu 
O E) t touj en 
10 es oujours 1 rieur = : 
Q, : ES TER 
kP 


m 


sur le cercle trigonométrique, les arcs et kA (au facteur 27 


Q @ 


près) diffèrent de moins de Loc . Lorsque m— ©, Q, aug- 
m m4] m A 
mente indéfiniment, au moins aussi vite que 2°, et Q = reste infé- 
À 1 m+l 
rieur à 2 * 
-l 
Donc ———— tend vers 0 et les points 
Que 
En 
ik 
Le \T m e2 iTkA 


sont très voisins. L'argument du terme de la suite er sl, tar 
+++ Qn_1) qui est le plus voisin de 1 est au moins égal à 


il Qn-1 27 Q n-1 
2 —_ — ————— = _ ———— 
(à Qu) 0 ( =) 


2 irmkA 


Finalement, le module de e - 1 est au moins égal à 
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à Q 
air m=-] 
(3 1-35 ) Q QT 
e ( A) on ed Re 0 ee 


Q Qns1/ Q Qui 


Comme Q 


m+l 


> 2Q,., cette quantité est elle-même supérieure à 


il 
. (8.12) 


C'est cette valeur que nous allons attribuer à la quantité n. 


Au total, l'erreur commise sur l'intégrale est majorée par 


Q 
2M ARTE 
+ e—— 
1Q7,  Nn 2 I 


Soit S la somme de la série |c,| +... + [ck|+... D'une façon 


générale, S est inférieure à c, ++ < MAT) es 1) 


n'est pas très utile, et peut en général être remplacée par une majora- 


, mais cette majoration 


2 
tion plus précise. Si par exemple À=1, ona |c,| us etS<MT. 


k?2 ? 6 
Quoi qu'il en soit, on obtient 
2M 45 
| 1, | < + nN 
ml 
c'est-à-dire 
8, 
Hé (8.13) 
avec 
2 M il 
œ_= ; n= 29Q,- 
Sd AuO Rs p . 


Nous allons interpréter la formule (8.13). 


À était positif, la suite a, décroît ettend vers 0. Pour une valeur 


re à . 1 
fixée de m, l'erreur |1,| décroît linéairement en fonction de N et tend, 


lorsque N—— w, vers une limite «, d'autant plus petite que m est 
plus grand. 


Etudions ce qui se passe lorsqu'on remplace m par m + 1. 


Dans le plan des (u, N), on remplace l'hyperbole H, qui a pour 
équation 
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8, 


UE 


par l'hyperbole H,. qui a pour équation 


m+l 


Or B,,, = 2SQ,,, est supérieur à $, et augmente indéfiniment avec m, 
pendant que a, décroît et tend vers 0. Les hyperboles H, et H,,,se 
coupent en un point qui a pour abscisse le nombre positif. 

be : B, 

XL À m1 


La disposition relative de H, et H,,, est celle qu'indique la 
figure 10 : 


Figure 10 
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Pour limiter l'erreur |1,|, on procède donc de la façon suivante : 
on choisit le nombre A. On calcule les Q, jusqu'à un rang m, tel que 
Œn, SOit suffisamment petit. On majore |I,| par 


Pm 
a, À FS 
Pour améliorer la précision, on construit, à partir d'un certain 
rang m,< m,, les hyperboles Ho, Hn,+1 ++, Hn.,. On conserve de cha- 
cune d'entre elles l'arc qui se trouve en-dessous de l'hyperbole précé- 
dente et au-dessus de l'hyperbole suivante. La courbe qui représente en 
fonction de N la limite supérieure de l'erreur commise est donc consti- 
tuée par une suite continue et décroissante d'arcs d'hyperboles, tendant 
vers 0 lorsque N —«, et qu'il est facile de construire si les Q, ont été 
calculés. Si l'on choisit pour A le nombre nr, qui n'est peut-être pas 
le nombre irrationnel donnant la convergence la plus rapide, on a 


Q m1 Q;, =7 Q, = 106 Q; = 113 
Q, = 33102 Q, = 33215 Q = 66317... 


[ll 


Y 


Il n'est pas sûr que cette méthode donne la meilleure évaluation 
possible de l'erreur commise. Elle est probablementassez grossière, 
car elle consiste à remplacer tous les dénominateurs e? ”*“- 1 par le 
plus petit d'entre eux, alors que certains ne sont pas très petits. Elle 
a l'avantage surles méthodes dontilsera question plus loin de ne partir 
d'aucune hypothèse sur le nombre irrationnel A. 


L'influence de la nature de l'irrationnel À peut s'interpréter de 
la façon suivante. Il y a intérêt à rester le moins longtemps possible 
sur une hyperbole H, fixée, et à passer à une hyperbole plus rappro- 
chée de l'axe des N. Il y a donc intérêt à ce que la suite B, soit régu- 
lièrement croissante, et cela de telle sorte qu'il faille changer le plus 
fréquemment possible d'hyperbole. Il semble que cette circonstance 
se produise lorsque À est un nombre algébrique, et plus particulière- 
ment un nombre quadratique. 


Si par exemple 


on a 
A =A +1, A =1+ 


Les quotients incomplets sont tous égaux à 1, et Q, est solution de l'é- 
quation aux différences finies 
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Q, = Qu, +0, 
avec 
Q, 4 ils Q, = 1 
On trouve 
m+l m+1l 
[9-0 os 
135 
Q; =2; Q, = 3, Qi=-39 4 2 ls ee 


m+]l . 
Lorsque m—>, É Æ tend vers 0. Ona donc approximativement 


m+l 


@ -Y5 À (8.15) 
Il existe deux constantes C et C' telles que, si m est suffisam- 
ment grand, on ait 


m 


C A 
IL ir + CN. 


(8.16) 


On vérifie que H, et H,,, se coupent en un point d'abscisse positive 
égale à 


- 1 Mm+1) 


(ex 
CG Te L 


>> 


dont l'abscisse augmente avec m. La disposition est celle de la fi- 
gure 11. 


IX - APPROXIMATION DANS LE CAS RATIONNEL - 


Si l'on remplace À par une fraction irréductible n ,; la suite des 


FE a : à 
réduites —- a un nombre fini de termes. A partir d'un certain rang, 


les termes ne varient plus. Tout ce qui a été dit au paragraphe 8 reste 
valable, sauf que la suite des hyperboles H, est limitée. Il y a une 
dernière asymptote qu'iln'est pas possible de dépasser. Cela ne prouve 
d'ailleurs pas que l'erreur commise ait une borne inférieure infran- 
chissable. Cela indique seulement que, par la méthode employée, on 
ne peut plus réduire indéfiniment cette erreur en donnant à N des va- 
leurs de plus en plus grandes. Ce dernier résultat est cependant exact, 
et nous allons maintenant expliquer pourquoi. 


L'erreur commise est égale à 
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Figure 11 


I SENS e À (9:47 


Si = est une bonne approximation de A, q est un nombre très 
grand, et N, bien qu'atteignant des valeurs très grandes, reste infé- 
rieur à q, et même petit devant q. Si le nombre k _ n estentier, cela 


provient de k, et non de n. Distinguons donc deux cas : 


1/ k est multiple de q. 


ik LB n N-1 
e "4 est alors égal à 1. La somme D) vaut 1. Le terme 


n=0 
correspondant de I, est exactement égal à c,. Dans le calcul de I,, on 
trouve donc la somme 
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CRTC TT Ce (9.2) 


qui est indépendante de N et irréductible. 


2/ kn'est pas multiple de q. 


Le coefficient de c, est égal à 


2 ik ÊN 
L Le (9.3) 
: imk Ë a 
Il est majoré par 
À L (9.4) 
2 Frrk 
CLR SE 


Il tend donc vers 0 comme è . On retrouveles résultats généraux 


de la discussion du paragraphe précédent. Mais on obtient ici une ex- 
pression de la limite de I, lorsque N——} w. On est sûr que la valeur 
absolue de l'erreur commise est au moins égale à 


ÉRIC Cr (9. 5) 
Si q est grand, cette quantité est naturellement très petite. 
X - RESULTATS DE RICHTMYER [4] - 


Richtmyer semble avoir été le premier à signaler, en 1952, 
l'intérêt éventuel de l'utilisation de la suite An dans les calculs d'inté- 
grales par la méthode de Monte-Carlo. Il a pu obtenir des majorations 
intéressantes de l'erreur I, lorsque À est un nombre algébrique. Voici 
quelques-uns de ses résultats (qu'il a réussi à généraliser aux inté- 
grales multiples). 


Rappelons d'abord deux théorèmes concernant les nombres 
algébriques. 


Le théorème de Liouville est relatif à l'approximation des nombres 
algébriques par les nombres rationnels. Supposons que A soit racine 
d'une équation algébrique à coefficients entiers de degré r (et ne soit 


pas racine d'une équation de degré moindre). Soit à une fraction ir- 
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réductible quelconque telle que |A - È 


<1, q > 0. Il existe alors une 


constante c, indépendante de p et q, telle que 
1 L 


C 
RRSNE (10.1) 


sauf peut-être pour un nombre fini de valeurs de p et q (l'inégalité est 
donc sûrement valable dès que q est suffisamment grand). 


Le polynôme dont À est racine se met en effet sous la forme 


BORA ON EE A CAN RE CE FAP H entier 
Donc 
e(B)-H(2-4 EE per 
q q q q 
Or 
JB-a,| -[B-a+A-A, < E-a|+ AR Ar A | 
q q q 
Les termes pe A; sont aonc majorés par des quantités indépen- 


dantes de : . Il existe un nombre c, dépendant seulement du polynôme 


e(®) 


g, tel que 


< £ je SA | 
c 1q 
D'autre part, 


Rap. entier | 
q q q 


L'entier du numérateur a un module au moins égal à 1, sauf si E est 
racine du polynôme g. Ce cas excepté, on a q 


») ee 
FOIE 
Finalement 
1 2 nee 
c1q q” 


D'où l'inégalité cherchée. 
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Dans le cas des nombres quadratiques, r est égal à 2, et le théo- 
rème de Liouville s'écrit 


| qA ele | (10.2) 


Sauf dans ce dernier cas, on peut beaucoup améliorer le théo- 
rème de Liouville. D'après le théorème de Roth, l'inégalité 


il 
|qA - pl TU (10.3) 
est valable quel que soit le nombre & > 0, aussi petit qu'on le désire, 


sauf peut-être pour un nombre fini de valeurs de p et q. On trouvera 
la démonstration de ce théorème dans le livre de Cassels [3]. 


Ceci posé, la méthode de Richtmyer consiste à partir de 
l'expression 


1 e2 iTKAN 1 
= Zu ———, (10. 4) 
k#0 e2TkA_ ] 


et à étudier d'une façon précise le comportement du dénominateur. 
On a 


les LES ain n KA: 


On remarque ensuite que : 


si O<x<z, Isinr x| >2x, 
son ; 
si 3 <x<1, IS rx 201 - <>) 


et que, d'une façon générale, 
(sin re 2 els (10.5) 
où l est l'entier le plus voisin de x. 
On a donc 
Le? TEE RAA 


D'après le théorème de Roth, on a, quel que soit 1, l'inégalité 
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1 


IkA -1| > 
kKô+1 


9 Kk#20; 


(sauf peut-être pour un nombre fini de valeurs de k). 


Il en résulte que, à partir d'une certaine valeur de |k|, 


2iTkAN 1 


e Ne 
EE —— 


2 


, (10.6) 


e2 iTKA = il 


et que la série dont I, est la somme est majorée par la série ayant 
pour terme général 


a. RE ER (10.7) 


Supposons pour fixer les idées que l'inégalité de Roth soit vala- 
ble sans restrictions. Dans le cas contraire, il faudrait isoler et cal- 
culer à part les termes, en nombre fini, qui ne sont pas justiciables 
de l'inégalité. On a alors 


1 


e 5+ 
ASS EuRe (10. 8) 
k=]1 


la somme étant étendue aux valeurs positives de k. 


Cette inégalité fournit une majoration de l'erreur I, lorsque la 
série du second membre est convergente. Or par hypothèse, 


M 
lex] +1 
Ik| 
d'où 
EC (10.9) 


Ily a convergence si À > 1 +6. Comm” 6 peut être choisi arbi- 
trairement petit, il suffitque À > 1. Le résultat de Richtmyer s'énonce 
ainsi de la façon suivante : 


THEOREME 10 - 


Si f(x) admet sur l'intervalle (0.1) un développement en série de 


Fourier dontles coefficients sont majorés par ceux de la série 1e où 
Ka 
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À estun nombre supérieur à 1, l'erreur commise sur le calcul de l'in- 
L 
tégrale f(x)dx en employant la suite arithmétique An (mod 1), où A 


0 ; 
est un nombre algébrique, tend vers 0 comme N Ce résultat est en 
RS ARS RP à | 


particulier valable si f(x), considérée comme une fonction périodique 
de période 1, admetune dérivée seconde continue, et une dérivée troi- 
= EE ———_— ———— — ——— — 
sième qui soit continue sauf peut-être en un nombre fini de sauts dans 


chaque période. 


Lorsque À est un nombre quadratique, on peut remplacer ô par 
0, et la démonstration de ce théorème ne fait pas appel au théorème 
de Roth, mais seulement au théorème élémentaire de Liouville. 


XI - QUELQUES RESULTATS SUR LE CALCUL ARITHMETIQUE 
DES INTEGRALES A L'AIDE DES POLYNOMES DE DEGRE 
SUPERIEUR AU PREMIER - 


Les méthodes utilisées pour étudier la suite An tombent en défaut 
dès qu'on remplace An par un polynôme de degré au moins égal à 2. 
On connaît assez peu de résultats précis sur le comportement de la 
différence 


TE 
[ ténéx - & E te) 


n=0 


lorsque x, est de la forme @(n) mod 1, où q est un polynôme ayant un 
coefficient irrationnel. Cette différence tend vers 0, et nous n'étu- 
dierons pas plus en détail comment elle tend vers 0. 

Nous nous contenterons d'examiner ce qui se passe pour 


p(n) = An? 


et lorsque À est remplacé par une fraction irréductible E. 
q 


Ilest commode de supposer que f(x) est développable en série de 
Fourier. Tout revient donc à l'étude du cas particulier où 


fee 0e k entier, 


et nous nous limiterons même à l'étude des sommes 


q 2 iTkB n2 


L 
€ à (11.1) 


n= 


Al 
o 


pour N = q. 
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Si k est multiple de q, ces sommes sont égales à 1, comme les 


= p 
Ne 7 précédemment étudiées. 


Si k n'est pas multiple de q, cen'est pas une restriction de choi- 
sir k = 1. Cela revient à faire jouer à la lettre p le rôle qui était pri- 
mitivement dévolu à kp. Considérons donc la moyenne 


1 ES 2 imBn? 
=, q 
te 


n=0 


S = (HIDE) 


Introduisons la fonction g(t), périodique, de période q, qui est 
égale à 


2 im Bn2 à 
e q SANTE NES 


C'est une fonction en escalier, développable en une série de Fourier. 
g(t) est de carré sommable, et la formule de Parseval lui est appli- 
cable. S est égale à la valeur moyenne de g{t) sur une période. Si donc 
ONÉCrIE 


Ce) 2im t 

g(t) = Ÿ ae * # (1123) 

on a 122 

Dia, 
On sait d'ailleurs que, si 

her iioss 

HSE 8 (t) g(t + hjdt (11.4) 
e 0 


désigne la ‘fonction de corrélation" de g(t), on peut représenter y(h) 
par une série de Fourier ayant comme coefficients les carrés des mo- 
dules de ceux du développement de g({t) : 


œ DETTES 


JE 
rh) 5 [affren (11.5) 
j=-0 
I1 en résulte que |a,|* est égal à la moyenne de y(h) sur une pé- 
riode : le carré du module de la somme $S est égal à la moyenne de la 


fonction Yy(h) sur l'intervalle 0 <h <q. 


On démontre d'autre part [1] que la fonction y (h) est continue et 
que, dans chaque intervalle (n, n +1), elle varie linéairement. Ses 
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valeurs pour h entier suffisent donc à la déterminer. Nous allons les 
calculer. 


Si h est entier, on a 


q=1 Zi [(n+h)22n2] 2m 2 q-1 Learn 
Ve 3 =e 1 NET (111261) 
pe di 


n=0 n=0 


ve 


© | 


L'expression de Y(h) est plus simple que celle de S. Il n'y figure ps 
qu'une progression géométrique. On trouve naturellement que y(h) = 
pour h=0eth = q. y(h) est nul pour presque toutes les autres ee 


de h. Il ne peut y avoir exception que si _ est entier. 


Si q est impair, cette circonstance n'a jamais lieu, et y(h) = 0 
pour tout h # 0 mod q. 


SI Ge on a 


im 1 2q°-1 in AL 
yta')=e oui DRE que : (112 
n=0 
y(q ) = 1 ou -1 suivant que q' est pair ou impair. 


Dans le premier cas (q impair), y(h) a pour moyenne à . 
Dans le second cas (q = 0 mod 4), y(h) a pour moyenne : : 


Dans le troisième cas (q = 2 mod 4), y({h) a pour moyenne 0. 
Il en résulte que, suivant les cas, |S| a pour valeur 
V2 
3 NC \a $ 
THEOREME 11 - 


air LB n2 


Le module de la somme S = ES 2 e 1" est égal à 


si q est impair 


si q = 0 mod 4 


siq = 2 mod 4. 


NOMBRES ALÉATOIRES, SUITES ARITHMÉTIQUES 


q impair 
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On remarque que la somme 


LR zirmBn 
DE 
nz0 


© | 


étudiée dans les paragraphes précédents, est toujours nulle. Ceci mon- 
tre que, lorsqu'onse limite aux sommes Ê S ,; l'approximation d'une 
intégrale à l'aide de suites Fi est moins bonne qu'à l'aide des suites 
Ph. Dans le premier cas, la différence entre l'intégrale et sa valeur 
approchée contient des termes qui sont petits avec + mais ne sont 


pas nuls, alors que ces termes disparaissent identiquement dans le 
second. 


Remarque. 
g=1 2i7 Ên2 F 
Les sommes Ye contiennent comme cas particulier, pour 


n=0 
p=1, les sommes de Gauss 


1 g-1 DATE n2 
— e q à 
L'État 
QT eZ 
e q 
En intégrant la fonction — sur un contour convenable du 
(OOMEE | 


plan complexe, on démontre que les sommes de Gauss valent 


ie 1 NUE 
ya ’ ÿq : » Va » 
suivant que q prend, mod 4, les valeurs 
0, 1 2e 3 


La méthode employée dans ce paragraphe ne fournit que les mo- 
dules des sommes de Gauss, qui suffisent d'ailleurs pour discuter le 
problème posé. 


XII - CONCLUSIONS - 


De l'étude qui précède, on ne peut pas tirer de conclusions ab- 
solues et définitives. Sous réserve de vérifications plus poussées, il 
semble qu'on puisse adopter les règles de conduite que voici : 


Pour calculer une intégrale par la méthode de Monte-Carlo, il 
y a intérêt à employer des suites arithmétiques plutôt que des suites 
de véritables nombres aléatoires. 
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Comme suites arithmétiques uniformément denses, on connaît 
celles qui résultent des théorèmes de H.Weyl et quelques suites de 
structure voisine. Les suites de polynômes p(n) mod 1, que nous avons 
étudiées, ne sont pas les seules qu'on sache construire, mais ce sont 
très probablement les plus maniables. Elles se prêtent bien au calcul 
électronique. 


Parmiles polÿnômes p(n), ceux qui donnent la meilleure appro- 
ximation sont les polynômes du premier degré An, où À est un nombre 
irrationnel. L'approximation est d'autant meilleure que la fonction à 
intégrer est plus régulière. On sait évaluer, assez grossièrement, 
une limite supérieure de l'erreur commise en remplaçant 


Qi TUST 
: f(x)dx par N 2 1(5), 
où x, = An mod 1, dès que f(x) admet un développement en série de 


Fourier dont les coefficients c, Sont majorés par 0e 
kK1+À 


Il semble qu'il y ait intérêt à choisir pour À non pas un nombre 
transcendant, mais un nombre algébrique, et même un nombre 
quadratique.(1) 


Si À est quadratique, et si À > 1, on sait que l'erreur commise 


est inférieure a où H est un nombre que l'on sait calculer. 


La méthode arithmétique s'applique au calcul des intégrales mul- 
tiples. On peut par exemple évaluer approximativement l'intégrale 


double 
[| f(x, y)dx dy 


ec 


étendue au carré 0 <x <1, 0 < y < 1, en la remplaçant par la somme 


le, 

N 2 1% Yn 
où la suite x,, y, est équirépartie dans le carré C. On sait qu'il en 
est ainsi si la suite 


Lx col 


(1) Par une étude plus détaillée de la structure du nombre irrationnel À, on peut 
sensiblement améliorer ces conclusions. Voir à ce sujet un prochain article 
de M. Bertrandias, où seront en outre développés plusieurs exemples numé - 


riques. 
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est équirépartie mod 1 pour tous les couples d'entiers I, 1, non simul- 
tanément nuls. 


Dans le cas de suites x, = An, y, = B,, cette condition est réa- 
lisée si les nombres irrationnels À et B ne sont pas liés par une rela- 
tion à coefficients entiers de la forme 


DASHIABI EI 


Exemple: 


A = V2, B = \3. 


La suite x, = An, y, = A(n +1), mod 1, ne satisfait pas à cette 
condition. 


Si p(n)est un polynôme de degré 2 au moins, ayant un coefficient 
irrationnel (autre que le terme constant), la suite 


x, ip (n), AN EEE À À mod 1 


est équirépartie dans le carré C. Elle peut servir au calcul des inté- 
grales doubles. 


La suite An semble la plus favorable au calcul des intégrales. 
Mais elle présente une corrélation interne. Si l'on a besoin de résou- 
dre par la méthode de Monte-Carlo des problèmes dans lesquels il est 
indispensable que la suite de nombres aléatoires ne présente pas de 
corrélation interne, il faut faire appel à des polynômes @(n) de degré 
2 au moins, par exemple An? mod 1. La précision des résultats est 
alors difficile à prévoir. On ne sait pas calculer d'avance le nombre N 
de termes qu'ilfaut calculer pour que l'erreur commise soit inférieure 
à une valeur imposée. 


La suite des parties entières mod 2 de p(n) est, du point de vue 
arithmétique, équivalente à une suite de nombres aléatoires prenant 
les valeurs 0 et 1, sans corrélation interne. 
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PROBABILITÉS À POSTERIORI DES PARAMÈTRES 
D'UNE LOI NORMALE 


Edgar BATICLE 


Je considère une loi normale dont la densité de probabilité est 
h/Vx e" fm?) On a effectué une épreuve qui a donné les résultats x: 
X2... Xn. Je Me propose de montrer qu'il est possible avec ces ré- 
sultats, supposés indépendants, de calculer la probabilité du couple 
aléatoire (M, H), correspondant aux valeurs certaines met h des pa- 
ramètres, sans qu'il soit nécessaire de faire appel au théorème de 
Bayes. 


Ce calcul repose sur la remarque suivante : ‘Si une variable z 
est une fonction croissante, ou décroissante, q(t) d'une autre variable 
tet siZ et T sont deux variables aléatoires telles qu'on ait Z = ®(T) 
on à : 


Pr{Z < z} Pr{p(T) < œ(t)} Pr Het 


ou 


JTE GE, 


PRES) 


Prip(T) < (t)} 
selon que la fonction (t) est croissante ou décroissante. 


Si, de plus la fonction (t)}, supposée croissante, varie de 0 
pour t = -œ à 1 pour t=+x, z = œ{(t) est la fonction de répartition 
de la variable aléatoire T. Si la fonction p(t) est décroissante et varie 
de +1 pour t = -w à O0 pour t = +, z = 1 - œ(t) est la fonction de 
répartition de T, et l'on a toujours Pr{Z < z}=7z:la distribution de 
z sur le segment (0, 1) est uniforme". 


Cela posé, je considère la fonction de répartition des n varia- 
bles aléatoires indépendantes X,, X,... X, c'est-à-dire l'intégrale 


multiple Z = ff... fer Er" (h/V 'dx, dx,... dx, dans laquelle les 
limites supérieures d'intégration sont x,, x,... x,. Soit F(x; x2 ... 
x,; M, h) cette fonction. Si les x,, x,... x, Sont donnés, on voit ai- 
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sément que F, est une fonction décroissante de met croissante de h. 


Je considère maintenant les variables-aléatoires Z, M, H cor- 
respondant respectivement à z, m, het telles que l'on ait pour les 
valeurs x, x, ... x, : Z = K(x, X, ... X M, H). 


1/ - Quelles que soient les valeurs attribuées à M et H la 
fonction F,(x,x,...; m, h) est une fonction de répartition des x et 
la coordonnée z est telle que PriZ < z}=2. 


2/ - Si on donne à H la valeur certaine h, on obtient 
Zh =, 20. 2 M, 7h) 
et l'on a pour la valeur certaine m : 
FC. XL here Pr 7 27; 
et, d'après la remarque préliminaire, on a : 
Pr{Z;<z}=/Pr(M > m | } 


De même, on aurait pour la valeur certaine m attribuée à 


nn Et) 


LAN EN(R Xe. x. 


m 
et pour H = h : 


HY(CexX 0. Xe Mon) zu pp 2,2 ep r IH che 

3/ - Mais, du fait que Z = F(x, x: ... xn, M, H}) établit pour 
toute valeur certaine z attribuée à Z une relation fonctionnelle 
entre M et H telle que si M croît, H décroît, et réciproquement, 


M > m entraîne H «< h. 


Par conséquent, la valeur commune de Pr{M > m |,} et 
Pr{H <h},} est égale à Pr{M > m; H <h}. 


Cela étant admis, on sait que si on appelle lu et 5 la moyenne 


et l'écart-type empirique des x, ... x,, l'intégrale multiple ci- 


dessus se ramène, à une constante près, à l'intégrale double 


anne 2, ee à » 
ffn'e nn loequ-n) 1 cr du do. Or, cette intégrale peut s'écrire : 


vu He tenir nn eee rats a(2) 
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Si, laissant fixes H et O, qui résultent de l'épreuve, je considère 
h et m comme variables, on pourra écrire : 


22 NN 
z = Kffn®e" BégemT n dh. 


Donc z - en supposant k déterminé par la condition que z = 1 pour 
m et h infinis (h étant d'ailleurs essentiellement positif) - est la 
probabilité que M est supérieur à m et H inférieur à h. 


La densité de probabilité du couple aléatoire (M, H}) est donc, 
à un facteur constant près : 


h°°! em *lot+çu-n)?] 


En intégrant par rapport à h de 0 à + on voit que la densité 


de probabilité de m est proportionnelle à [o? +(u - m)]”?, et si 
on pose be = cot$ on trouve que la densité de probabilité de 


8 est A,sin"#8, avec A, J sin?8 = 1. D'où A, =1/YT.T, :T.. 


2 2 


De même, on trouverait en intégrant par rapport à m que la 


2 


densité de probabilité de H est K e"" f?hr-2 avec K 7 rc ous 


Il est à remarquer que l'expression trouvée pour la densité 
de probabilité du couple (M, H}) est celle qu'on aurait obtenue en 
appliquant le théorème de Bayes avec la probabilité a priori vir- 


tuelle ee proposée par divers auteurs. 


Les résultats précédents trouvent une application particulière- 
ment importante dans l'interprétation des résultats d'épreuves 
quantitatives sur échantillons, quand on admet la validité de la loi 
de Laplace -Gauss. 


1/ - On peut, en premier lieu, se proposer de rechercher 
quelle est la probabilité a posteriori d'une valeur donnée de la 
moyenne (m) du caractère étudié dans la population-mère. 


La formule énoncée ci-dessus donne le résultat cherché. On a : 
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ue TG) ps 
> n 
EciMen) er $ d®, 
re 
formule dans laquelle n est l'effectif de l'échantillon et $ un angle 
(nm -u) 


auxiliaire tel que l'on ait cot’#$ = TNT 


C'est d'ailleurs le problème inverse qu'on a le plus souvent 
à résoudre : on cherche la valeur de m qui a une probabilité donnée 
de ne pas être atteinte ou d'être dépassée. 


I1 s'agit là d'une technique connue sur laquelle il n'y a pas 
lieu d'insister et qui rejoint la théorie de l'intervalle de confiance. 
D'aiileurs cot $ n'est autre que la variable t de Student. 


2/ - I1 faut remarquer que la probabilité de la moyenne n'est 
pas toujours la plus importante. Par exemple, quand on se préoc- 
cupe du degré de sécurité sur lequel on peut compter dans l'utili- 
sation du lot d'où l'on a extrait l'échantillon soumis à l'épreuve, 
ce qui importe c'est que le caractère envisagé ait une probabilité 
d'être atteint (ou, selon le cas, dépassé) qu'on puisse considérer 
comme négligeable. 

Pour répondre à cette question, nous cherchons la densité de 
probabilité d'une valeur x de la variable, compte tenu des résultats 
de n essais. 


Cette densité de probabilité est, d'après ce qui est établi ci- 
dessus, proportionnelle à l'intégrale 


PfR EN OR re ARR 
prise de 0 à +w pour h et de -x à +« pour m. 


Posons 


et (net) rn =D x, + x 
Le coefficient de - h° dans l'exponentielle devient alors : 
CAE TR RENE T PR PR 1e 


Faisons le changement de variable 
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= (n+ re [(m - u')? Te TU er 
L'intégrale double devient : 
© +0 NARUE Dinant 
LEE ein 1) L)0+ sie. L'AlL2Zadu dr 


On voit qu'elle est proportionnelle à 


n + 
ft - u'Ÿ Free u'*] Zäm 
-@ 


qui devient si on fait le nouveau changement de variable 
m - u' = t(s'? "x HUANIGR 
rs nt £ n 
FREE AREAS Re 
Er ) 


Elle est donc proportionnelle à (s'? - ae 2, 


E valuons s'? - u'? en fonction de o,u, x. On a 


SUD our - nl {(n +1Èx? + x°) - [2x,)* + 2x 2x, + x) 
Er [(n + 1{no° + nu° + x?) - (n?u° + 2n x + x°)] 


nn Lin Lot (x Lie 


2 

n 2 (x -H) 

(nel) |s+En 
Gus. 


Si on pose rot = cot‘p, la densité de probabilité de 


sera À, sin" et sa fonction de répartition 


n 
['a sin" do = — Jet” 32 


La loi de probabilité de ® est donc ,1a même que celle de $ 
et l'on passe de $ à q@ par la relation cot 79 = Ga 4 1) cot” p. Une 
seule table suffit pour le calcul des probabilités de m et de x. 


Si on appelle k (K = cot}) l'écart normé de x par rapport à 
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la moyenne empirique p, les valeurs de k, en fonction de n pour 
les probabilités 0,001; 0,002; 0,010, sont données par le tableau 


ci-dessous : 


1042" 15547 


3 4 5 6 7 8 9 


SOINS NES 
4,56 4,453,9 3,6 ... 3,09 
51€ 547 40/40 


M.Fdgar Baticle m'a remis cette note le 7 mat 1960, à la fin d'une 
séance de Séminaire de Statistique. Bien que tous les participants du 
Séminaire n'aient pas été d'accord avec les conclusions de l'auteur nous 
publions ce travail en raison de son intérêt. 


Daniel Dugué 
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